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Abstract 

This paper is a continuation of [Brlvr] and [Ivrl, Ivr5, Ivr6, Ivr7]. I derive 
sharp spectral asymptotics (with the remainder estimate 0(h l ~ d + / u r /i 1 ~ r ~ l? ) for 
c?-dimensional Schrodinger operator with a strong magnetic field; here h and fi are 
Plank and binding constants respectively and magnetic intensity matrix has constant 
but not full rank 2r at each point and q = d — 2r. 

Introduction 

0.0 Preface 

In this paper I consider multidimensional Schrodinger operator 
(0.1) A = A + V(x), A =J2 Pj9 jk (x)P k , 

j,k<d 

Pj = hDj - fjVj(x), h e (0, 1], fi>l. 
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It is characterized by magnetic field intensity matrix (Fjk) with Fjk = (d x .Vk—d Xk Vj), which 
is skew-symmetric d x d- matrix, and (F£) = (gi k )(Fk p ) which is unitarily equivalent to skew- 
symmetric matrix (g^ k ) z{Fjk)(gi k ) 2 . Then all the eigenvalues of {F ] k ) (with multiplicities) 
are ±i/ p (/ p > 0, p = 1, . . . , r) and of multiplicity q = d — 2r where 2r = rank(i^). 

In this paper I assume that magnetic field intensity matrix has a constant but not full 
rank: 

(0.2) rank(F iJfe ) = 2r = d - q, q > 1, \ ({F^)^)' 1 ] < Co 

where Af = Ker(Fjk) and Af ± = Ran(Fjfc) and (under certain conditions) I derive sharp 
spectral asymptotics (with the remainder estimate 0(h l ~ d +n r 'h l ~~ r ~ q ). The typical example 
(already studied in [Ivrl, Ivr6]) is 3-dimensional magnetic Schrodinger operator. 
Condition (0.2) is equivalent to 

(0.3) rank(F jfc ) = 2r = d - q, \f p \ > e p = 1, . . . , r. 

I consider operator in some domain or on some manifold X with some boundary con- 
ditions, assuming that it is self-adjoint in L 2 and I denote by e(x, y, r) Schwartz' kernel of 
it's spectral projector. 

There is a profound difference between full-rank case considered in [Ivr7] and non-full- 
rank case considered here: in its canonical form full-rank magnetic Schrodinger operator is 
reduced to the family of r-dimensional fi~ l h-pdos and therefore we can expect asymptotics 
with the remainder estimate (in the smooth case) as good as C(fih~ 1 ) r ~ 1 x (l + (/i/i)~ r ) x 
( K H~ 1 h 1 ~ d + /i r ~ 1 /i 1 ~ r ) where the first factor is the remainder estimate for a single operator 
yU -1 /?,-pdo and the second factor is the number of these operators to be taken into account; 
however one needs some non-degeneracy (microhyperbolicity) condition; otherwise the re- 
mainder estimate could be as bad as C( K jih 1 ~ d + ji r h~ r ^j which is the magnitude of the 
principal part as /i > hr x . 

On the other hand, in its canonical form non-full-rank magnetic Schrodinger operator is 
reduced to the family of g-dimensional Schrodinger operators which are also r-dimensional 
/i^/i-pdos and in the best case we can expect asymptotics with the remainder estimate 
Ch x ~ q x (fih" 1 ) 1 " x (l + (/i/i)~ r ) with the principal part of magnitude h~ q x (iih~ 1 ) r x (l + 
(fih)~ r ); however, non- degeneracy condition is not needed at all (if smoothness allows) as 
q > 2; case q = 1 is more complicated. 

We see that there is a more subtle distinction between q — 1 and q > 2; actually q = 2 
might stay apart from q > 3 due to the lack of the smoothness. There is even more subtle 
distinction between r = 1 and r > 2 but it makes the difference only as q = 0, 1 setting 
d — 2, 3 a bit apart (and in fact it plays role mainly in the non-smooth case). 
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0.1 Assumptions. I 

Note first that M = Ker F is integrable foliation of q- dimensional spaces, which means that 
for any y there exists g-dimensional manifold N = N y 3 y such that T X N = Af(x) at any 
point x G N. Then there exists a coordinate system x = (x"'\ ar 1 ) = (x±, . . . ,x g ; x q+ ±, . . . , x&) 
with x = const along leafs N y of the above foliation. This (almost) fixes x and also im- 
plies that in this coordinate system Fjk = as j = 1, . . . , q (or k = 1, . . . , q). Then there 
exists function such that Vj = djcj) for j = 1, . . . , q and one can eliminate Vj for all 
j = 1, . . . , q by a gauge transformation. So, we can assume that Vj — as j — 1, . . . , q and 
then Vj = Vji^x- 1 -) as j — q + 1, . . . , d. 

Now we assume that these reductions are already made and thus 

(0.4) F jk = 0, Vj = 0, djV k = Vj = 2r + 1, . . . , d Vfc 

and after these reductions 

(0.5) gi k ,Fj k EC l '°, VeC 1 '*, 

(0.5)' q = 1 =>► gi d = S jd , 

(0.6) = djcf>j, ^ E C 1 ^ 

where (7, a) >z (I, cr) b (1, 1)- We also assume that 
(0.7) \V\ > e 

Strange-looking condition (0.6) will be needed in some cases to make global reductions to 
a canonical form. 



0.2 Canonical Form 

It is known (see e.g. [Ivrl], Chapter 6) that if X = M d and gi k ,Fj k ^V are constant then 
the magnetic Schrodinger operator is unitarily equivalent to operator 

(0.8) f S (h 2 D* + £ h*ti} + V 

l<j<r r+l<j<d-2r 

and that e(x, x, r) = £^ /IW (r) where 

u d ^ r {2-K)- d+r ^h- d+r Y, (r-J2(2a J + l)f jf ih-vfj l ---f r ^ 

a£Z+ r j 



(0.9) £r(r 



def 
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and Uk is a volume of unit ball, g = det(<7 jfe ) _1 . 

Now in the general case (for variable V and, may be, g jfc , Fj k ) I am interested in the 
asymptotics of J e(x, x, 0)ip(x) dx as h — > +0, // — > +00 where ^ is a fixed function, smooth 
and compactly supported in X. 

It is known [Ivrl] that in the smooth case for d = 3 one can reduce magnetic Schrodinger 
operator to its canonical form 

(0.10) /i 2 a mtP j tk (x2,x 3 ,fi~ 1 hD 2 )x 

m+p+j+k>l 

x (x 2 + t T 2 h 2 I%) m {hD z )* l r^*- 2 ^ k h k 

with the "main part" 

(0.10)o ai,o,o(«2, a? 3 , n~ 1 hD 2 ) (fx 2 xj + h 2 D\) + /i 2 L> 2 + a ,i, (^2, a; 3 , fx~ l hD 2 ). 

with a^oo = F o \J/ ) 00,1,0 = ^ ^0 an d certain diffeomorphism \l/ : ^l^ 1 x K — > IR 3 ; 
we ignore terms with fc>lorm + p + j>2. Then one can replace harmonic oscillator 
(h 2 Df + n 2 xf) by one of its eigenvalues (2a + l)fih. 

Multidimensional case is much more tricky. The main part of the canonical form will 

be 

(0.11) Z l a Wi x '"> fi~ 1 hD")Z j + 

j,k<r 

h 2 DjbJ k (x", x'", yT x hD")D k + a%(x", x'", fi^hD") 

2r+l<j,k<d 

with Zj = (hDj — ifiXj), Z* = (hDj + ifiXj) j = l,...,r, x" = (x r+ i, . . . ,x 2r ), x'" = 
(2r + 1, . . . , d). Here (a,jk)j,k=i,...,r is a Hermitean positive definite matrix with eigenvalues 
fx, . . . , f r and (bjk)j,k=2r+\,...,d is a re& l symmetric positive definite matrix. We would like 
to have (aj k ) in the diagonal form; then instead of (0.11) we would have 

(0.12) ^2fJ r (a/ , y,n- 1 hD")(H j -nh) + 
3 

^DjbJ^x", x"\ ii~ l hD")D k + a (x", x"\ [T x hD"Y 

2r+l<j,k<d 

with 1-dimensional harmonic oscillators Hj = (h 2 D 2 + fi 2 x 2 ). However second-order res- 
onances fj = f k with j k prevent us; in the general case we can assume only that 
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locally 

(0.13) £ £ a, fc (x",x" / ,/i- 1 /i J D / ')Z;Z fc + 

^ h 2 DjbJ k (x", x'", fi^hD")D k + ao(z", fi^hD", x") 

2r+l<j,k<d 

where m G DJt are disjoint subsets of {1, . . . , r} and eigenvalues of each of matrices (&jk)j.kem 
are close to one another. This leads to the necessity of the matrix rather than the scalar 
non- degeneracy (microhyperbolicity) condition. 

The reduction of non-principal terms of operator are prevented by higher-order reso- 
nances (see [Ivr7]); we discuss it later. 

0.3 Microhyperbolicity 

Thus we arrive to the final microhyperbolicity condition at point z = (x\x , f) and thus 
at point x = (z): 

Condition 0.1. Assume that 

(0.14) First of all we assume that set {1, . . . ,r} is split into disjoint groups m G fJUl(z) 
and matrix (a^) is block-diagonal and in a small vicinity of z its blocks are close to scalar 
matrices f m I# m where / m are disjoint. 

(0.15) Finally, we assume that for each real vector r = (r m ) me sr»t such that 
| J2 m em Tm + ^1 — e there exists vector £ = £(z, t) G M. 2r+g such that 

(0.16) a 2 Y,(*( a ik% 1 ))C]a > ei VC G C 

as long as 

(0.17) a m =J2 ^fcCjCfc = tfr m Vm G 971. 

To impose this condition we need to assume first that symbol a does not vanish; thus 
we need assumption (0.7). 

Note, that contrast to [Ivr8] I no longer consider 91 partition of {1, . . . , r} joining in one 
group all possible third-order resonances as well; the reason is that now we need to consider 
evolution with \t\ < eo instead of \t\ < eofi as it was in [Ivr8] and therefore third-order 
resonances are not a problem anymore here. 
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Condition 0.2. The above microhyperbolicity condition should be modified as \ih > €q: 
in (0.16) — (0.17) one needs to consider only r m = ^j Gm (2o!j + l)/i/i/i with a» G Z + ; now r 
runs a discrete set. 

On the other hand, for constant matrix (Ft) condition (0.16)-(0.17) is equivalent to 
(0.18) |W|>ei; 
and for fih > eo one should modify it to 

(0.19) |5^(2ai + l)/A/i + V| + |W|>ei Va G Z +r . 

i 

Then we basically do not need condition (0.7). Also in this case one can take £ which does 
not depend on r. 

0.4 Results 

The first group of results uses no canonical form and is related to a relatively weak magnetic 
field. We consider two cases: without and with microhyperbolicity condition respectively. 

0.4.1 Weak Magnetic Field. I 

Theorem 0.3. Let conditions (0. 1) - (0.5), (0.7) be fulfilled, fx < h s ~ l and 
(0.20) Either q > 2 and (I, a) >z (1, 1) or q = 1 and (/, a) >z (1, 2). 

Then there exist two framing approximations 1 such that for both of them 

(0.21) U = \ J^e(x,x,T)-S MW (x,T)^(x)dx\<Ch 1 - d + Cfih 1 - d (fih\\ogh\)i. 

Remark 0.4. (i) Theorem 0.3 uses no non-degeneracy condition and it does not benefit 
from stronger smoothness assumptions. 

q 

(ii) For ji < jii(q) = C(h\ log h\)~^ remainder estimate O(h ) holds; in particular, 

11 3 

p>i(i) = C(h\ log/i|)"S, /i 1{2 ) = C(h\ log/i|)~2 and fi 1[g) > C(h\ \ogh\)~s as q > 3. 
1 I.e. operators such that A + > A > A~ in the operator sense. 
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0.4.2 Weak Magnetic Field. II 

To push \x in theorem 0.3 up 2 we need microhyperbolicity condition and extra smoothness. 

Theorem 0.5. Let conditions (0.1) - (0.5), (0.7), (0.20) be fulfilled and fj, < h 8 " 1 . 
Let microhyperbolicity condition 0. 1 be fulfilled. Then 

(i) The following estimate holds 

(0.22) K < Ch l - d + Ch~ d (fih) l+ i | log h\ l+ %- a + Cfih- d (fih) 1+1+ i | log h\ l+1+ 2 ~° ; 

(ii) Furthermore, if either [i < fa — C(h\ log/i|)~^ or £± = in microhyperbolicity 
condition, then one can skip the last term in the right-hand expression of (0.22): 

(0.23) H < Ch 1 -' 1 + Ch- d (fih) l+ i | log h\ l+ i- a - 

(Hi) In particular, 1Z < Ch 1-11 as q = 1, (I, a) = (|, 1) and \i < fa- 
Note that for fi close to h^ 1 the last term in (0.22) could be larger than the principal 

part. 

0.4.3 Intermediate Magnetic Field. I 

We assume first that microhyperbolicity condition 0.1 holds. 
Let us assume now that 

(0.24) fa = Ch~*\ log/i| - 3 <^<fa = eh~ l \ \ogh\~ 1 . 

The low bound is justified as q > 2 since otherwise theorem 0.3 provides 0(/i 1_d ) remainder 
estimate. For q — 1 justification is that we will need to take e = Cfa 1 and it is smaller 
than e = Cfa\\ log/i| used in the previous results if and only if \i > fa. The only exception 
(when we can take e = C(fa 1 h\ log /i | ) 5 ) is the case r = 1 but then q > 2 since we assume 
that d > 4. 

Theorem 0.6. Let conditions (0.1) — (0.7), (0.20) and (0.24) be fulfilled and let microhy- 
perbolicity condition 0. 1 be fulfilled. Then 

(0.25) K < Ch l ~ d + Ch~ d (fih\ \ogh\)2fa l \ \ogh\~ a . 

Corollary 0.7. Under microhyperbolicity condition for \i < fa sharp remainder estimate 
K < Ch^ d holds as either q> 2, (I, a) h (1,1) or q = I, (I, a) h (§, 1). 

2 Actually we need to do this for q = 1 only. 
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0.4.4 Intermediate Magnetic Field case. II 

Now let us consider intermediate magnetic field without microhyperbolicity condition. Then 
essentially we have a degenerate matrix operator. To such operator we can apply rescaling 
technique in its crudest form leading us to 

Theorem 0.8. Let conditions (0.1)- {0.7), (0.20) and (0.24) be fulfilled. Then 

(i) For q > 3, (/, a) = (1,1), (Z, cr) = (2, 1) estimate TZ < Ch 1 - d holds; 

(ii) For q = 2, (/, a) = (1,1), (I, cr) = (2, 1) estimate 

(0.26) TZ < Ch l - d + Cfihl~ d 

holds; in particular, TZ < Ch 1 ~ d as ji < h~% ; 
(Hi) For q = 1, (I, cr) = (1,2) estimate 

(0.27) TZ < Ch l ~ d + Cfihl~ d + Cfih 1 -' 1 

holds. 

Note that estimates (0.26), (0.27) are better than (0.21) with q — 2,1 respectively iff 
/i > h~z \ \ogh\~2 . 

The second terms in (0.26), (0.27) are due to remainder estimates for (/-dimensional 
Schrodinger operator (it is not important that operator in question is also r-dimensional 
/i^/i-pdo). These estimates correspond to the C 1 ' 1 -smoothness and cannot be improved in 
the general case even if we assume much larger smoothness because we have a matrix, not 
a scalar h-pdo as for d = 3. However, in certain cases we can do better than this; we need 
to consider q — 1,2 only. 

In (0.27) the third term estimates an approximation error arising when we skip irre- 
ducible 0(fi~ 1 ) terms due; we can prevent these terms under certain conditions. 

Theorem 0.9. Let q = 1, (1,2) ^ (I, a) H (2,0) and conditions (0.1) - (0.7), (0.24) be 
fulfilled. Then 
(i) Estimate 

(0.28) TZ! = \ J(e(x,x,T)-S MW (x,T)-S^(x,T)^(x)dx\ 

< Ch^ d + Cfihl~ d + C^-h 1 -^ log/i|-f 

holds where 

(0.29) \S^(x,t)\ <C^h 1 ~ d 
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will be defined by formula (5.67) and in the general case one cannot improve estimates 
(0.27), (0.29); 

(ii) If one of conditions 

(0.30) There are no third-order resonances 

(0.31) gi k = const, Fj^ = const 

is fulfilled then estimate (0.28) holds for 1Z. 

As r > 2 condition (0.24) is fulfilled let us consider 
(0.32) n (x, h, t) = #{a e Z +r , ^(2^- + l)fj(x)h + V(x) < r}. 

3 

Let us assume that estimate 

(0.33) |/ (n (x,h,r') - n (x,H,r')) dx\ < Cfr r (\r-r'\ +v{h)) 

JB(0,1) 

\/hE (0,1) Vr,r : |r| < e, |r'| < e 

holds with u(h) = o(H). 

Estimate (0.33) holds with u(h) = h for sure and it holds with u(h) = h K , k > 1 if either 
fj have constant multiplicities and rank{V^-, . . . , Vj^} > k — 1 for all x or = const 

and system . . . , ^} has some Diophantine properties 3 . Under condition (0.33) we can 
improve theorems 0.8, 0.9 as q = 1, 2: 

Theorem 0.10. Let q = 1,2 and conditions (0.1) - (0.7), (0.24), {0.33) be fulfilled. 
Then 

(i) For q = 2, (/, a) = (1, 1), (/, a) = (2, 1) estimate 
(0.34) ft < Ch l - d + Cu(fih)h 2 " d 
holds; 

(ii) Forq = 1, (1,2) ^ (/, tr) r< (2,0) estimates 

(0.35) ft < C7i 1 - d + C7/-3/i3-rf|iog/ l |§ +Cv(fih)(h 1 ~ d + n~h~ d ^, 

(0.36) ft/ < C/i 1 ^ + C^h~ l hh- d \ log h| 

C(y{nh) x ^3- d + /i-5/ l - d |log/ i |-f)/ i - d 

3 See examples 4.13, 4.14 for details. 
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hold; 

(Hi) For q = 1, (1,2) ^ (I, a) ^ (2,0) as one of conditions (0.30), (0.31) is fulfilled the 
following estimate holds: 

(0.37) n ^Ch 1 -' 1 + Cn*- l h?- d \logh\*-' 7 + 

Cu{fih) (h l ^ d + /j-h- d \ log/i|-i) . 

Improved version of theorem 0.10(i) when for TZj estimate (0.37) holds will be formulated 
later (theorem 5.12). 

0.4.5 Intermediate Magnetic Field case. Ill 

Assume now that 

(0.38) Matrix (F£) has eigenvalues ±if m of constant multiplicities. 
Then microhyperbolicity condition 0.1 is equivalent to 

(0.39) \V+Y1 frt\ + MV+ fi T i)\>£ Vn, . . . ,r r e E + . 

Further, if /2//1 = const, . . . , f r / f\ = const then (0.39) is equivalent to 

(0.40) |V^| > e. 

Ji 

Theorem 0.11. Let conditions (0.1) - (0.1), (0.20) and (0.38) be fulfilled. Then 

(I) For q = 2, (1, 1) ^ (I, a) ^ (2, 0), (I, a) = (2, 1) estimate 

(0.41) 1l<Ch l - d + Ciih l - d+ &\\ogii\-^) 

holds. 

In particular, 1Z < Ch l ~ d as (I, a) = (2,0). 

(II) As q = 1, (1, 2) ^ (I, a) ^ (2, 0), (I, a) = (2, 1) estimate 

(0.42) ft/ < CTi 1 "'' + C>/i 2 ~ d ~^| log/ipl^ + Ch 1 -*^ log/i|) V~'l log/i|- CT + 

Cy-^ 1 -^ log/i|-f + C/i^/it-^l log/i|-s 

holds. 



INTRODUCTION 



11 



(in) Asq = l, (1,2) r< r< (2,0), (Z,tr) = (2,1) and one of conditions (0.30), (0.31) 
is fulfilled estimate 

(0.43) ft < C , /i 1 " d + C , /i/i 2 ~ d ~^|log/i|"(^ + 

Ch 1 ~ d (fih\ log/i|) V -, | log/il - " + C^-h 1 -*] ]ogh\S 

holds. 

Note that the right-hand expression of (0.42) contains an extra term (in comparison 
with (0.43). 

Also note that for q = 1 estimate (0.41) is better than (0.21) as 
(0.44) // +1 /i|log/i| CT+1 > 1; 

in particular, as (/, a) = (2, 0) we get estimate 1Z < Ch l ~ d + C ' [ihh~ d which always is better 
than (0.21) as /i > 

Under condition (0.33) one can improve theorem 0.11: 

Theorem 0.12. Let conditions (0.1) - (0.7), (0.20) and (0.38), (0.33) be fulfilled. Then 
(i) Asq = 2, (1, 1) ^ (I, a) r< (2, 0), (I, a) = (2, 1) estimate 

(0.45) TZ < Ch^ d + Cu(fj,h)h 1 - d+ '^\ \ogh\~^ 

holds. 

(n) Asq = l, (1,2) ^ (I, a) ^ (2,0), (I, a) = (2,1) and one of conditions (0.30), (0.31) 
is fulfilled estimate 

(0.46) K <Ch^ d + Ch- d (fih\ \ogh\)^- l \ log/i|- CT 

Cu(fih) (h 1 ~ d -^\ + Cfi-h- d \ \ogh\~^ 

holds. 

Improved version of theorem O.ll(ii) will be formulated later (theorem ??). 

0.4.6 Strong and Superstrong Magnetic Field 

Let us assume that 

(0.47) e/z" 1 ! log/il" 1 < n < ChT 1 . 

Then we must pick up 

(0.48) e > C(^~ l h\ log/i|)a + Ch\ \ogh\ x h\ \ogh\ 

Still theorem 0.8 remains true: 
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Theorem 0.13. Let conditions (0.1) - (0.7), (0.20) and (0.47) be fulfilled. Then 

(i) For q > 3, (I, a) = (1,1), (I, o~) = (2, 1) sharp remainder estimate 1Z < Ch x ~ d holds. 

(ii) For q = 2, (I, a) = (1, 1), (I, a) = (2, 1) estimate (0.26) holds: 71 < Cfihl~ d . 
(m) Forq = 1, (1,2) ^ (I, a) ± (I, a) = (2, 1) estimate (0.27) holds: 71 < Cfihl~ d . 

Theorems 0.6, 0.9 and Corollary 0.7 require some kind of localization with respect to 
a m defined by (0.17); however, this localization is not needed if either r = 1 and d = 4, 
or microhyperbolicity condition (namely condition 0.1 as fi < eh" 1 or condition 0.2 as 
eh" 1 < fi < Ch" 1 ) holds with vector i which does not depend on r. If this is the case, 
theorem 0.6 and corollary 0.7 still hold. 

This localization with respect to be done on Hermite decomposition level as long 

as there no higher order terms; we can skip them with error not exceeding 0(/i 1_d ) as q = 2. 
Then for q = 2 theorems 0.6, 0.9 and corollary 0.7 remain true. 

Also these statements remain true for q — 1 if there are no unremovable 0(fi~ 1 ) terms 
i.e. if one of conditions (0.30), (0.31) is fulfilled. 

If q — 1 and there are unremovable 0(/i -1 ) terms, this procedure brings an approxima- 
tion error 0(n~^h^~ d ) under microhyperbolicity condition and 0([i2h 1 ~ d ) without it. This 
error estimate would not exceed remainder estimate as fi > eh" 1 . 

Further, in this particular case when q = 1, r > 2, /i< efih and both conditions 
(0.30), (0.31) are violated we can still localize a m by microlocal arguments similar to those 
of [Ivr8] provided 

(0.49) e > C( f i- 1 h\logh\)K c ^ loshl + Ch\\ogh\. 

One can check easily that an approximation error is less than the corresponding remainder 
estimate unless I = 1 and we are looking at theorem 0.6; in the latter case an approximation 
error is 

(0.50) C(fi~ 1 h\logh\)h Cllhlloehl h" d \\ogh\- a . 
Thus we arrive to 

Theorem 0.14. Theorem 0.9 remains true with condition (0.24) replaced by (0.47). In 

particular, in estimate (0.28) the third term is less than the second one. 

Theorem 0.15. Let conditions (0.1) - (0.7), (0.20) and (0.47) be fulfilled, q = 1,2. Let 
microhyperbolicity condition 0.1 (as fih < e) or 0.2 (as \ih > e) be fulfilled. 

(i) Let either q — 2, or ji> eh" 1 or microhyperbolicity condition be fulfilled with vector 
i independent on r , or one of conditions (0.30), (0.31) be fulfilled, or I > 1. 
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Then estimate 7Z < Ch l d holds, 
(ii) Otherwise estimate 

(0.51) 71 < Ch^e ^ ^ log/il 1 -" 

holds. 

The same " drop unremovable terms" arguments could be applied in frames of theorems 
0.10-0.12 leading us to 

Theorem 0.16. Under condition (0.24) replaced by (0.47) 

(i) Theorems 0.10, 0.11(i),(iii), 0.12 remain true; 

(ii) In frames of theorem 0.11(H) estimate 

(0.52) 71 < Ch x - d + Cfih 2 - d -Th\ log/i|~^W + Cph 1 ' 4 

holds and one can skip the last term unless (l,c) >- (2, —2). 

0.4.7 Ultrastrong Magnetic Field 

Finally let us assume that 

(0.53) pL > Ch~ x 

Then we need to modify the second part of condition (0.5), replacing (0.5) by 
(0.5)* g>'\ F jk e W = V + (Ah fj e C l >° 

3 

and assuming that 

(0.54) W + 2f j /ih>e Vj. 

Under these assumptions we have really a scalar operator, and should take 
(0.55) e = Ch\\ogh\ 

and microhyperbolicity condition 0.2 transforms into 
(0.56) \W\ + |W| > e. 
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Theorem 0.17. Let conditions (0.1) - (04), (0.5)*, (0.6), (0.7), (0.20), (0.53) and (0.54) 
be fulfilled. 
Then 

(i) The principal part of asymptotics is 
(0.57) J £f w (x,T)tf}(x)dx = 

UJ d -2r(27T)- d+r ^h- d+r f(r-J2 f0h -Vfh--- f r ^(x) dx 

of magnitude (fih) r h~ d . 

(ii) For q > 3, (I, cr) = (1,1), (I, cr) = (2, 1) estimate 

(0.58) K < C(fih) r h 1 - d 

holds. 

(ii) For q = 1,2, under microhyperbolicity condition 0.2 (coinciding with (0.56)) esti- 
mate (0.58) holds. 

(m) For q = 2, (1, 1) ^ [I, a) ■< (2, 0), (7, a) = (2, 1) estimate 

(0.59) TZ < C(fih) r h 1 ~ d + C(fih) r h 1 - d+ ^ | log h\~%* 

holds. In particular, for (Z, cr) = (2,0) estimate (0.58) holds, 
(w) For q = l, (1, 2) ^ (/, a) < (2, 0), (I, a) = (2, 1) estimate 

(0.60) K < C(^h) r h}- d -m\\ogh\~^m 

holds. 

0.5 Plan of the paper 

Tiny section 1 contains preliminary remarks. In section 2 we develop weak magnetic field 
approach to cover both weak magnetic field case and the exterior zone in the intermediate 
magnetic field case. In section 3 we reduce operator to its canonical form. Section 4 is 
devoted to remainder estimate with implicitly defined main part of asymptotics, which is 
simplified and calculated in section 5. Tiny section 6 contains some closing remarks. 
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1 Preliminary remarks 

We repeat section 1 of [Ivr6]. First of all, proposition 1.1 of [Ivr6] remains true (arguments 
repeat those of subsection 1.1 of that paper). On the other hand, proposition 1.2 is replaced 
by two following statements with 



(1.1) £ mw (t) = 



(2n)- d+r ^h- d+r V / A ■ • ■ frV$dt 

-^^J{T-EA2a j +l)f jf ih-V+\e»\>><0} 



where fj are eigenvalues of matrix (F^) = (g^ k )(F km ) and gi k , F km , V are corresponding 
mollifications; here e depends on 

Proposition 1.1. Let q > 2 (I, a) = (1, 1). Let (I, a) = (2, 1) in statements 
(i) Let p, < h~2 . Let us define 



;i.2) 



{Cp 1 h\\ogh\ p > pi = max(/i 1 , (ph\ \ogh\) 2 . 

Cph\\ogh\ 



P < Pi 



with p = Then the following estimate 

(1.3) TZ'=\ J ^(x)(^ MW (x,r)-£" MW (x,r))dx| < Ch 1 ' 4 V|r| < e. 

(ii) Let hr*\ \ogh\~ 2 < p < h" 1 ] \ogh\~ 1 . Let us define 



;i-4) 



Cp 1 /i|log/i| + ^ ViP s P > Pi = log 
C/i" 1 p< (nh\logh\)* 



with p = Then estimate (1.3) holds. 

(Hi) Let n > h~ l \ \ogh\~ 1 , e = (p,~ l h\ \ogh\)* . Then 

(1.5) W < C(l + ph) r h l ~ d V|r|<e. 

Proposition 1.2. (%) Let q = 1, (1,2) -< (I, a) -< (2,0). Let (I, a) = (2,1) in statements 

(a), (m). 

(i) Let fx < h~3. Let us define e by (1.4). Then estimate (1.3) holds. 
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(ii) Let h i | log /i | 3/i < h 1 \\ogh\ 1 . Lei us define e by (1-4) • Then under microhy- 
perbolicity condition 0.1 

(1.6) n' < Ch 1 -' 1 + Ch~ d {fih\ log h\)^~ l | log h\~ a V|r|<e 
and without it 

(1.7) ft' < C/i^ + C/i-^/z/illog/iD^-'llog/il-^ + C^-^^llog/il-f 

V|r| < e. 

(mj Lei /i > h" 1 ] \ogh\~ 1 , e = (n~ l h\ \ogh\)^ . Then under microhyperbolicity condition 

0.1 

(1.8) W <C(l + iih) r h l ~ d V|r|<e. 
and without it 

(1.9) ft'< C(l + /i/i) r ^" d +(l + /i/i) r rV^og/»r f V|r|<e. 

2 Weak magnetic field case 

Analysis in this case follows one section 5 of [Ivr6] without serious modifications. We assume 
that ji < h 6 ^ 1 . Let us assume that f operator conditions (0.2)-(0.5) are fulfilled. 

2.1 Preliminary remarks 

Let us start from the case (further restrictions to follow) 

(2.1) e>C7i|log/i|, (i<C- 1 h- 1 \logh\- 1 . 
Obviously 

(2.2) [P j ,P k ]=i f ihF jk 
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and all the commutators 

(2.3) h-^Xj] = J2(9 jkp k + Pk9 jk ), 

k 

(2.4) h~ l [A, P m ] = h' 1 Pj[9 Jk , P m ]P k + h~ 1 [V,P m ] m = 1, . . . , 2r, 

(2.5) fi-'h-'lA, P m ] = Y,{9 jk F jm P k + P J g jk F km ) + 

^ l h- 1 PA9 3k , P m ]P k + ^h-'iV, P m ] m = 2r + l,...,d 

jk 

have symbols of the class T l,a in any domain {a(x, £) < c/i -2 } of the bounded energy where 

(2.6) a(x,£) = a (x,0 + V, a (x, f) = /i 2 ^g^p^x, £)p k (x, £), 

Classes JF Z,<J are introduced in [Ivr6]. To finish preliminary remarks let us notice that the 
following statement holds 

Proposition 2.1. 4 Lei / E C£°(M. d ), f = 1 in 5(0, Co) wit/j /arye enough constant Co. 
Let T > h 1 " 6 . Then under condition (2.1) with large enough constant C = C s 

(2.7) \Ft^ h -iT{l-f(pi(x,Z),P2(x,Z),...,Pd(x,0)) X T (t)u\<Ch s 

Vx,yeB(0,^) Vr<c. 

One can see easily that symbol f(pi(x, (,),p 2 (x, £), . . . ,Pd(%, £)) ^ s quantizable under 
condition (2.1) and the proof of proposition 2.1 is rather obvious. 

2.2 Analysis in the exterior zone 
Proposition 2.2. 5 Let p < eh' 1 ] \ogh\~ 1 6 and let 

(2.8) M> sup 2V/( t + £ 
4 Microlocal boundaries for Pj. Compare with proposition 5.1 of [Ivr6]. 

5 Finite speed of propagation with respect to x; cf. Proposition 2.1 [Ivr6] and Proposition 1.1 [Ivr7]. 
6 Lower bound condition is not needed in this statement as well as non-degeneracy conditions. 
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with arbitrarily small constant e > 0. 

Let (f>i be supported in B(0, 1), <p2 = 1 in B(0,2), \ be supported in [—1, 1] 7 . Finally, 
let f = Ch\ \ogh\ <T <T[ = e . 

Then 

(2.9) \F t ^ h -i T x T (t)(l -<fa,MT(x-x))</> ljMT (y-x)u(x,y,t)\ <Ch s 

Vx, y Vr < ei 

where here and below e\ > is a small enough constant. 
Proposition 2.3. 8 Let /i < eh 6-1 and let 

(2.10) f = Che' 1 ] \ogh\ <T<T[ = e . 
Then 

(2.11) |F t ^-i T XT(*)(l-02,MT(^ 

< <7/i s Vx, y Vr < ei 

where here and below ip is an admissible function supported in B(0, |) and ip x = ip(x), 

Proposition 2.4. 9 Let p < eh 5 - 1 and condition (2.10) be fulfilled. Let 

(2.12) p > C/i£ _1 |log/i| + Cp~\ 
T/ien /or 

(2.13) C(/i _1 /i| log/i|)5 + Cp- 2 /i| log/i| <T<ep 
inequality 

(2.14) |F t ^-i T XTWrg^w(x, <C7i s Vr< ei 

/io/ds /or Q y = <p lM T(hD'>' + £'")• 

7 I remind that we pick up all such auxiliary functions to be admissible in the sense of [Brlvr] . 
8 Finite speed of propagation with respect to P'"\ cf. Proposition 5.2 [Ivr6] . 
9 Singularities leave diagonal; cf. proof of Proposition 5.3 [Ivr6] . 
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Now let us pick up 

(2.15) e = Cp~ 1 h\\ogh\. 
Then any T G [e/x _1 ,ep] satisfies condition (2.13) as long as 

(2.16) p> pi = Cfi^ 1 + C{fih\ log h\)^. 

This is a definition of the exterior zone and e is defined by (2.15) there. 

Then the contribution of the exterior zone to an approximation error does not exceed 

(2.17) C /(^^)'(l + ^^)|log/ i |-V- 1 ^x^ 

J P P 

with all integrals here and until the end of the section taken from pi to 1, and (2.17) is 
0{h l ~ d ) as either q > 2, (I, a) >z (1, 1) or q = 1, (£, <r) >: (1,2) (as stated in propositions 
l.l(i),1.2(i)). 

On the other hand, Tauberian arguments exactly as in [Ivr6] imply that the semiclassical 
remainder estimate is 

(2.18) Ch^J ^jpi-Up 

or Ch^ d f p q ~ 2 dp which is 0{h}- d ) as q > 2. 
So, we have proven 

Proposition 2.5. As q > 2, e = Ch\ \ogh\p~ 1 , (I, a) = (1,1) contributions of the exterior 
zone > pi} with p > pi = CpT x + C(/ih\ log/i|)2 to &o£/i approximation error and 

semiclassical remainder estimates are 0(h 1 ~ d ). 



2.3 Exterior zone: special analysis for q = 1 

As g = 1 (2.18) with T{ = ep results in Ch l ~ d \ \ogp\ in and to get rid of | log/i| we need 
increase T[. We consider few cases in the increased generality and complexity. 

(i) If gi k and Fjk are constant it is easy: we can take T[ = Cp\ logp| 2 as in section 5 
of [Ivr6] provided (l,o~) >z (1,2) and we can consider propagation in an appropriate time 
direction (in which ^ increases); see section 5 of [Ivr6] for details. Namely, let us cover 
B(0, 1) by balls of radii 7 = p Sl (with small exponent 5i > 0). In each element one of the 
following cases holds: 

(2.19) \d Xd V(x,0\ x g > g = Collogpl" 2 , 

(2.20) \d Xd V(x,£)\<Q. 
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During propagation the classical trajectory does not leave the initial ball. Furthermore, in 
case (2.20) |^| = |9 X(J V| does not exceed g and |^| > \p — g\t\ > ~p as \t\ < T[. 

On the other hand, in case (2.19) |^| is larger than g and does not change sign for 
\t\ < T( . Then |^| >\p + g\t\ for < ±£ < T[ with an appropriate sign. 

In both cases shifts away from its initial position and for time t (of the appropriate 
sign in the case (2.19)) this shift is of magnitude p\t\ + g\t\ 2 as ep~ l = T < ±t < T[ (with 
an appropriate sign) and we can take a prescribed T[ = ep\ logp| 2 and recover 0(h 1 ~ d ) in 
(2.18). 

(ii) Similar arguments work if fj have constant multiplicities. Really, we can rewrite 
then symbol as 

(2-21) a {x^)=e d + p 2 Yl Wv-i+vD + V 

with rjk being linear combinations of pi, ■ ■ ■ ,P2r (with JF 1,2 coefficients) such that 

(2.22) p{r] k ,r] m } = A km 

modulo linear combinations of r]i,...,r]2r with T 0,2 coefficients where A^ m = 1 ii k = 
2j — 1, m — 2j for some j, Ak m = — 1 if k = 2j, m — 2j — 1 for some j, Ak m = otherwise. 
Then derivative of p 2 (r]2j-i + iV2j)(V2k-i ~ ^2fc) along classical trajectories 

(2.23) \p 2 {a, (7 ?2i _ 1 + iV2j)(V2k-i ~ iV2k)}\ < M(l + p\f, - f k \) Vj, k 

is bounded as fj = f k and during time T its variation does not exceed CT which is way 
less than eg. This could be justified on the quantum level as well because one can quantize 

symbol / ^f?" 1 ((?72j-i + iT)2j)(V2k-i ~ ^?2fc) _ A) j with function / supported in [— c, c] and 

|A| < c. 

Consider now {a, modulo O(p) it is equal to — d Xd V + p 2 J2j k a jk r lj r ]k with JF 0,2 
coefficients. Then we can find a', k such that 

(2.24) {a, £ d - n ^ a jkVjVk} = ~d Xd V + p 2 ^ a "k(V2j-i + ir}2j)(r}2k-i ~ ir)2k) 

j,k j,k:fj=fk 

modulo 0(g): we remove all non-resonant terms but previously we make //-mollification 
of Oij k thus making 0(g) error. But then all the above arguments work since correction 
Z 1 J2j u a jkVj r lk does not exceed Cp~ x < ep. The role of — d Xd V is played now by right-hand 
expression of (2.24). 
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(iii) Consider now the general case. First of all let us make again 7 = p Sl -covering and 
consider one element. Then the symbol could be rewritten in the form 

(2.25) a(x, f ) = i 2 d + ^ fl m + V, a m = p 2 ^ a jk{V2j-i ~ iV2j){V2k-i + ir]2k) 

m j,k£m 

with Hermitean matrices a m = (cijk)j,kem whose eigenvalues are fj(x), fj and differ by less 
than e (more than 2e) as j, k belong to the same group m (different groups, respectively). 

Consider traceless matrices b m = a m — tr a m ■ I m with corresponding unit matrices 
I m and matrices c m = d Xl b m calculated at some point of this element. Without any loss 
of the generality one can assume that each c m is a diagonal matrix and with eigenvalues 
classified into groups: Xj and differ less than eg (more than 2eg) as j, k belong to the 
same subgroup n C m (different subgroups, respectively). Then one can cover each element 
of partition by 71 = p 5l+&2 subelements of two kinds: 

(a) The total measure of this type elements (of all elements together) does not exceed Cg 
and thus their total contribution to the remainder estimate does not exceed Ch l ~ d \ logp|~ 2 
which after integration over y gives 0(h 1 ~ d ). 

(b) With fj, fk different by more than 71 as j, k belonging to different subgroups. Then 
we can rewrite (2.25) with n (which indicates finer partition) instead of m) and to prove 
that I {a, a n }\ < C^~ &i where £3 could be made arbitrarily small (by taking all the previous 
exponents small). Then a n are conserved modulo 0(g) in the classical evolution. 

Moreover, we can make even finer subpartition depending on subelement such that 
\fj — fk\ < Cg if i, j belong to the same element of it. 

But then f n fi 2 X^en^ij-i + vij) are conserved modulo 0(g) and then it is also true for 
a n and a m as well and also 

{a, £ d ~ V a JkVjVk} = Y A j(^2j-i + vlj) mod 0(g). 

j,k j 

But then arguments of (i),(ii) work and the total contribution of elements of type (b) does 
not exceed Ch l ~ d \ logp|~ 2 as well and after integration over -£ we get 0(h 1 ~ d ) again. So 

Proposition 2.6. As q — 1, e — Ch\ \ogh\p~ 1 , (1,<j) = (1,2) contributions of the exterior 
zone > pi) with p > pi = Cp~ x + C(ph\\ogh\)^ to both approximation error and 

semiclassical remainder estimates are C^/i 1 ^). 

Similar arguments were applied in [Ivr5] to study propagation near the boundary. 
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2.4 Main theorems. I 

On the other hand, the standard estimate in the complimentary zone — Pi} brings 

remainder estimate Ch}~ d x pp\ x Ch l ~ d x p(ph\ log ^ | ) 2 + Ch 1 ~ d [i 1 ~ q . 
Therefore we have 

Proposition 2.7. Let conditions (0.4) — {0.7) be fulfilled an let either q > 2, (/, a) >z (1, 1) 
orq = l, (I, a) >z (1,2). T/ien 
(%) For < e/i 5-1 estimate 

(2.26) ^ = |r(e^)(r) - /T 1 / F Hfl -VXr(t)r(^)| < 

C7i 1_d + C(i(jJth\ log/i|)2/i 1_d Vr : |r| < e x 

holds with T = T = efi" 1 . 

(ii) Furthermore, this estimate holds with any T e [T ,T] with T = Ch\ \ogh\. 

Statement (ii) follows from the fact that in the standard settings F t ^ h -i T XT(t)F (uip) is 
negligible for indicated T and r (I remind that \ is supported in [— 1, — ~] U [~, 1]) and 

(2.27) h' 1 I F t ^ T , X T(t)T(mP) dr' = iT~ l F t ^ h -, T x T {t)T{uij) 

J —00 

with x(t)=t~ 1 x(t). 

Proposition 2.7 implies Theorem 0.3 after we prove (by successive approximations or 
reference to standard results rescaled) that 

(2.28) h- 1 I F t ^ T ,x T (t)T(u^)dT 

equals J S MW (x,T)ip(x) dx modulo remainder estimate which is easy. I leave it to the 
reader. 

On the other hand, this approach is valuable even if theorem 0.3 is not our target 
because \i is pretty large. Namely, from the analysis of two previous subsections follows 

Proposition 2.8. Let conditions (0.4) — (0.7) be fulfilled. Then 

(i) for p, q < ji < eh' 1 ] \ogh\~ 1 with small enough constant e estimate 



(2.29) K 1Q = \r(Qei>)(r) - h- 1 f F t ^ h -x T ,Qx T {t)T{u^)\ < 



Ch l ~ d q>2, 
Ch l - d \\ogh\ q = l 
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holds as |t| < e± where Q = I—f(hD"') with admissible function f supported in < 2pi} 
and equal 1 in < pi} with pi = C(p~ 1 h\ log/i|)2 ; T = e\pT l . 
(ii) Further, ifq=l, (I, a) >z (1j2) and either 

(a) n < h 5 - 1 , or 

(b) fj have constant multiplicities, or 

(c) microhyperbolicity condition 0.1 is fulfilled on level 
then TZiq < Ch l ~ d as |r| < e\. 

We will use this proposition to handle exterior zone in the case of the intermediate 
magnetic field. 

2.5 Main theorems. II 

Further, applying technique of subsection 7.2 of [Ivr6] and of proposition 1.4 of [Ivr7] we 
can prove easily 

Proposition 2.9. Let p, q < /j, < eh~ x \ \ogh\~ 1 and microhyperbolicity condition 0.1 be 
fulfilled. Then for Ce~ l h\ \ogh\ <T<T[ = e\ estimate (2.14) holds. 

Selecting in the exterior zone {p > p{\ 

(2.30) e = Cfih\ log/i| (^-Y + Cp^h\ log/i| 

P 

and combining with proposition 2.4 we get then that for Q with the symbol supported in 
the exterior zone estimate (2.14) holds for e/i" 1 < T < €\. 

In the complimentary zone {p < p{\ we can select e = Cp>h\ \ogh\ (as /i > h~z \ \ogh\~*) 
and then for Q with the symbol supported there estimate (2.14) also holds for e/i^ 1 < T < 
T[ = €\. Thus we arrive to 

Proposition 2.10. Let h~z \ \ogh\^^ < p < eh~ l \ \ogh\~ 1 and microhyperbolicity condition 
0.1 be fulfilled. Let e be defined by (2.30). Then estimate (2.14) holds with Q = I and any 
Te[f ,T{],f = e^- 1 ,T{ = e 1 . " 

Combining with the standard results rescaled and applying standard Tauberian tech- 
nique we arrive to estimate (2.26) again. 

On the other hand, an approximation error does not exceed 

(2.31) C(nh\logh\) l+ i\logh\- a h- d 
as long ass < C/i -1 (i.e. for /i < Ch~^\ logh\~^). 
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Further, for /i > Ch 2 1 \ogh\ 2 an approximation error does not exceed 

(2.32) Cp(ph\ log h\) l+l+ i I log h\-~ a h~ d + C{ph\ log /i|)' + 2 1 log h^'hr*. 

However if £± = we can take mollification with e = Chp~ l | log h\ with respect to all 
variables x and also an extra mollification with e defined by (2.30) with respect to x'"; then 
mollification error will be given still by (2.31). 

This implies Theorem 0.5 after we prove (by successive approximations) that expression 
(2.29) equals J £ MW (x, 0)ip(x) dx modulo remainder estimate. This easy part of the proof 
I leave to the reader. 

3 Canonical form 

3.1 Intermediate magnetic field case 
3.1.1 Approximation 

We can assume now that 

(3.1) fti {q) = C{h\ \ogh\Y^ <p< eh- 1 ] \ogh\~ 1 

(a case of larger p we'll consider later) and from this section point of view there are two 
different cases q — 1 and q > 2; under microhyperbolicity condition q — 1 falls in the latter 
as well. Anyway we need to consider zone Z x = {|n < pi = C(fih\\ogh\)2} and in this 
zone in the general case we already set e = C{p~ l h\ \ogh\)^. However we reset it to a far 
larger value e = Cfi' 1 thus defining 

£ _ (Cfi- 1 ^) 3 + Cp- 1 h\\ogh\ asp>p!, 

y CpT x as p < pi 

with a possible increase later. 

For q > 2, p, > /x 2 = h~^ \ log h\~^ (and thus pr l \ log h\~ a < ph even as (/, a) = (1, 1)) 
this will bring an approximation error not exceeding 

Cp r h- r (ph\\ogh\) q 2- 1 p- l \hgh\- a x ( /i /i) 1 - r |log/i| 

where the first factor estimates an error in each term in the sum defining £ MW and the 
second factor estimates the number of the terms effected; the result is 0{h l ~ d } for sure. 
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As q = 1 let assume that 10 
(3-3) p 2 d ^ (»- l \loghD* < Pt=W- 

Consider strip {p < < 2p} with p 2 < p < pi where an approximation error can be 
estimated similarly by 

Chr d pp\ + Cph x - d p- x p\ 

and after summation over p we get that the contribution of of zone {p2 < < p±} to an 
error does not exceed 

(3.4) Ch~ d (ph\ \ogh\)^p~ l \ \ogh\- a + Ch}- d p}- l *\ log/i|-f . 

In the zone {|£ w | < C/O2} no more than C(ph) l ~ r terms in the sum defining £ MW are 
affected and an error in each does not exceed Cpfh~ r p2\ this brings an approximation error 
estimate Cph 1 ~ d p2 which is the exactly the second term in (3.4). 

Note that (3.4) is O^ 1 ^) for all p > Chr\ \ log h\ 4 iff (I, a) h (2,1). 

On the other hand, if q = 1 and microhyperbolicity condition 0.1 is fulfilled an the 
approximation error does not exceed 

(3.5) C/i- d (/i/i|log/i|)^- i |log/i|- CT 

and it should be weighted against theorem 0.5; so in this case for p < P2 = Ch~^ \ \ogh\"^ 
we use theorem 0.5 and for /i > ^ we will use the intermediate magnetic field construction; 
an approximation error estimate in both cases is 0(h l ~ <v ) for all corresponding p iff (/, a) y 




Proposition 3.1. Let e be defined by (3.2). Then 

(i) As q > 2, (Z,cr) y (1, 1) and p > h~z \ log h\~2 a mollification error is 0(h 1 ~ d ); 

(ii) As q = l, (l,o-) y (1,2), p > h~~i I log/i|~3 and (3.3) is fulfilled a mollification error 
does not exceed (3.4); 

(ii) As q = 1, (I, a) b (1 ; 2), p > h~z \ \ogh\~2 and microhyperbolicity condition 0.1 is 
fulfilled a mollification error does not exceed (3.5). 

10 This is not necessarily true as q = 1; however without this assumption no improvement of the previous 
section results is possible. 
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3.1.2 Reduction. Main part 

So we start from operator 

(3.6) A= Pj9 ik Pk + V : Pj^hDj-^Vj, 

l<j,k<d 

(3.7) [Pj, P fc ] = ifihF jk , 

(3.8) Fjk = as j > 2r + 1 or k > 2r + 1; rank(i^-fe) :7i fc = i j ... ) 2 r = 2r, 

and we can assume without any loss of the generality at some point x matrix (Fj^) is in 
the canonical form; namely 

{fj as j = 1, 3, . . . , 2r - 1, fc = j + 1, 
as j = 2,4, . . . ,2r, fc = jf-l. 

otherwise. 

Without any loss of the generality we can assume that in its vicinity 

(3.9) F jk = ^ j < 2r, k < 2r and 3m e 971 : j, fc e m. 
Then following [Ivr6, Ivr7] we can prove 

Proposition 3.2. Let 

(3.10) C = fi 2 (J2P Jk PjPk) W 

j,k<d 

with ft k = f3 k i G T l '° ' . Then modulo operators with symbols belonging to n~ 2 J rl ~ 1 '' J 

(3.11) e^- lh - Hc ( Pj ye-^ lh - ltC = {^{e M )]p kj 

k 

with matrix A = {h.jk)j,k = (Fj P )((3 pk ). 
One can see easily that due to (0.4) 

e ={o I 

with A' = (Aj fc )j fc < 2r , A" = (Aj k ) j<2r,k>2r+i- It provides us with the first step of the 
reduction: 
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Step 1. Applying transformation exp(|/x 1 h l C) with L defined by (3.10) with appro- 
priate coefficients (3^ k we can transform operator A to the block form 

(3.12) V 2 J2^ WajkW ^ + ^ E Pja ikw Pk+<% + \{B + B*), 

j,k<r j,k>2r-\-l 

(3. 13) B = (V V km PjPkPm + VPS) 

j,k,m j 

(3.14) Q = pj + ipj+r] a jk ^0 => 3m : j, k G m 

where I remind that DJl is the resonance partition at point x. 

Here and below a jk G -T 7 ^, a G ^' CT , # H G .F^ -1 ' ", 6* G ^- 1 - <7 and while a J ' fc = a jk (x) 
we have ^ fe ' = & kl (x, /i" 1 ^); ^ = b^(x, /i" 1 ^) and c w means h- rather than /i~ ^-quantization 
of symbol c. 

Step 2. Applying an appropriate gauge transformation we can achieve 

(3.15) = unless k = r + 1, . . . , 2r and = Vk(xi, . . . , x 2r ) for = r + 1, . . . , 2r. 

where condition (0.6) provides C' +1 '°-regularity of new V^. This is the only reason for this 
condition. Then 

(3.16) Pj = Vj=r + l,...,2r r, A ; (./•".(') Vj = l,...,r 

with x' = (xi,...,x r ), f = (6, ■■■,&), a/' = (x r+ i, . . . ,Z 2 r), f" = (£r+l, • • • , 6r), Aj G 

Following [Ivr6, Ivr7] and applying transformation e* M £ with £' = Yli<j<r Aj-Dj+r we 
can transform operator (3.12) to the form (3.12)-(3.13) but now with 

(3.17) Q = Vj + iVj+r, Vj = Vj+r = Xj j = l,...,r 

Vj = j j = 2r + 1, . . . , d. 

However in contrast to the previous step of the reduction now 

a? k = a jk (x", p-'C; x"'), a = a (x", fT 1 ?'; x'") 
rather than a? k = a^ k (x), ao = ao(x). 

Remark 3.3. (i) One can see easily that if fj have constant multiplicities, operator can be 
reduced on Step 1 to the form 

(3.18) » 2 Y,$ W f7<7 + V 2 E Pja jkw Pk+< + \{B + B*) 

j<r j',fc>2r+l 
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and this form will be preserved after Step 2. 

(ii) One can see easily that if q = 1 then after Step 2 by means of transformation e lh c 
with C" = ^(hDdk w + k w hDd), k = k(x", Xd, Ai _1 £") one can reduce operator to the same 
form (3.12) but with a dd = 1. 

(iii) In the general case we can achieve a? k = 5jk for j, k — 2r + 1, . . . , d only assuming 
that gi k , are constant. 

3.1.3 Reduction. Junior terms 

I remind that junior terms are of the form \{B + B*) with B defined by (3.13) with 

Remark 3.4. As (1,1) di (^, cr) -< (2, 0) we can rewrite B in the same form but with 

(3.19) V km = V km (x", iT 1 ?'; x m ), V = V(x", aT 1 ^; d") 

modulo operator B' with the symbol belonging to fi~ l \ log/i| _CT jF ' . Our reduction is done 
modulo operators of this type. 

Then applying operator e lfJl 2fl l£l with 

(3.20) A = ( J2 ^"W)*, 

a:l<|a|<3 

7] a = rji 1 ■ ■ ■rj^ d , one can reduce operator to the same form (3.12) with 

(3.21) £=//(E V km Q}Um.+ E V km Vmt]Qk+ 

j,k,m<r j,k<r, m>2r+l 

E ^' fcm ^w) w + ( E 6 ^) w 

j,k,m>2r+l j>2r+l 

Here in the first term V km = unless \ fj — fk — f m \ < e and in the second term V km = 
unless \fj — fk\ < e (so there are the third and second order resonances respectively). 
Further, V km = always for constant g^ k , Fj^. 

Also, in (3.20) (3^ belongs to C^ 1 '" 7 , C l ~^ a as | a | > 2, \a\ = 1 respectively; also 

Jjjkm g (jl—l,cr ^ £jl—l,tr 

Remark 3.5. (i) In the intermediate magnetic field case |?7j| < C/i _1 in the microlocal 
sense (on the energy levels below c) and therefore junior terms are of magnitude and 
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skipping them we would arrive to an approximation error Ch~ d (/j,h\ \ogh\) 2 (as either 
q > 2 or q = 1 and microhyperbolicity condition holds) which is 0{h 1 ~ d ) as q > 3 and 
0(h l ~ d \ log/i|) as q = 2. So, in the case of an intermediate magnetic field only for q = 1,2 
we will need to take junior terms into account at all. This will be a minor annoyance 
for q = 2 but a major obstacle for q — 1, when without microhyperbolicity condition this 
error would beO(V d 0/i|log/i|) 2 ji 1 + ji 2 h 1 d ) but even with the microhyperbolicity these 
terms cannot be simply ignored. 

(ii) Actually, not all of these terms are equally bad: terms containing factor(s) r\ m with 
2r + 1 < m < d are of magnitude 0(fi~ 1 \C,"'\) and can be simply ignored as q = 2 but even 
as q = 1 they are minor annoyances. Only first term in (3.21) are really important. 

(iii) One can see easily that if q — 1 and fj have constant multiplicities, then by means 
of transformation e iM h Cl with 

m jr',fcGm 

with c jfc = c> k (x", Xdi /i" 1 ^"), Cj = Cj(x", Xrf, /i" 1 ^") one can remove from B all terms con- 
taining factor r] d thus arriving to 

(3.22) B = ^ 2 ( ^ m CjCfeCm) W 

So, we constructed operator T with the properties described below: 

Proposition 3.6. Let (l,cr) = (2, 1), (1, 1) ^ (l,o~) ^ (2,0). T/ien inere exists an operator 
T such that 



(3.23) T*TQi> = Qijj, 

(3.24) T*Q{T{I - Q)$ = 0, 

(3.25) T*(I -Q 2 )TQ^j = 0, 

(3.26) T*ATQ^ ~ AQ^, mod log/il""^ ' ) 



as ^ G Co°(-B(0, 1)), <Q = Q(hD"') has a symbol which is supported in < 2(7^} 

and equal to 1 in {|£ w | < 2Cp{\, Qi and Q 2 are the same type operators but with symbols 
which are supported in < 2pi} and < 2C 2 pi} ^ nd e Q ual to 1 in < pi} 

and < C 2 pi} respectively with large enough constant C, A is the reduced operator 

(3.12), (3.21). 

As q = 1 and fj have constant multiplicities, operator A is given by (3.12), (3.22). 
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Remark 3.7. Note that while (3.23)-(3.25) hold modulo negligible operator, (3.26) is sat- 
isfied modulo /i~ 1 /i-pdos with symbols of pT | log h\~ a T°'°, and to emphasize the difference 
I use "~" instead of "=". 

In contrast to almost all paper in propositions 3.6, 3.8 b w means yU _1 /i-quantization of 
symbol b. 

We also need the following properties of operator constructed above: 

Proposition 3.8. Furthermore the constructed operator has the following properties: 
For functions $ = 1 m B(0, f) and $ E C£°(B(0, \) 

(3.27) i)Qi) = T*^ w TQij, 

(3.28) i) ~ VoO", x'", C) + Yl ^( x "> mod loghl^F ' 

l<j<d 

and for operator Q = Q w (hD"') as above 

(3.29) T*QT^ (g(/(x,0))^ mod/i-'llog/^r^T, 

(3.30) f(x,0-C- Yl «i^-^ _1 «o mod/i-^log/ir^ ' 

l<j<2r 

with coefficients Kj = Kj(x", x'", p^ 1 ^") and all estimates holding for \r)j\ < Cp~ l , ipj,ctj G 

Proof. We prove more difficult (3.29), (3.30). 

(i) On Step 1 of the main part reduction we considered operator e^ lh ltc with the 
corresponding symplectomorphism which is equivalent to the map f t : (x, •••,%) —> 
(x, rji, . . . , rjti) which is of T l ' a class. Let Q t = (Q(ft)) j then the principal symbol of 
operator 

(3.31) ^ft + ^lAft] 

vanishes and then as long as (7, cr) b (2,1) and £ > C(p~ 1 h\ log/i|)i the norm of this 
operator does not exceed 

Cfi^h- 1 x (-) 2 < C^/if/z/illog/i)- 1 = C^llog/il" 1 
P 

for p > C(ph\ \ogh\)^. Then the norm of operator 

(3.32) ^(e-^-HLQ^h-H^ 
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does not exceed Cp 2 | log/i| 1 and then the same is true for operators 
(3.33) e-^ lh ' Hc Q t e^ lh ' Hc -Q and e^ h ' Hc Qe-^ h ' Hc - Q t . 

(ii) The same arguments work perfectly on Step 2 of the main part reduction. 

(iii) Now consider the junior terms reduction. First of all, there are removable terms 
J2i<j<2r fitfj w ith (3j = Pj(x", p" 1 ^"), j3j G T l ~ 1,a which are removed by a transformation 
of the same type as above with £ = p~ l Xa<j<2r P'jVj an d P'j = p 1 ^,"), G T l ~ x ' u as 
well. 

Then the corresponding Hamiltonian map is of class JF /,<J and then the principal symbol 
of operator (3.28) vanishes and the norm of this operator does not exceed 

C/i^llog/ir x/i- 1 /^ 1 x ((-) 2 + p^h(-) x^-'llog/iH <Cfi- l \\ogh\-° x ^ 

V v p v p J p z 

which is less than Cp~ l | log h\~ u and then the same estimate holds for operators (3.32)- 

(3.33) . 

With the removable third and second order terms the construction is the same but 
simpler since these terms are coming from the principal part which belongs to T 2 ' 1 for 
sure. □ 

Remark 3.9. In [Ivr6] in the case d = 3 we tooke = C(p~ 1 h\ log/i|)5 instead of Cp~ l . The 
above construction would not work. In this case note first that the rescaling x h- > e~ 1 {x — x) 
and taking h = p~ x he~ 2 = \ogh\~ l we will see that the original operator x{hDs/p) 
becomes x(hD 3 ) and p~ x h~ l C becomes vt{x, hD) with v = p~ 2 e~ l = C^ 1 (p 3 h\ \ogh\)~2 < 
ei and therefore h c becomes h-pdo with "analytic" symbol (more precisely its symbol 
and those of Q satisfy [Brlvr]). Then e~ lfl h c Qe % ^ h c becomes ft-pdo which means that 
in the original scale it is /i" 1 /i-pdo with JF°'° symbol supported in {|£ 3 | < 2Cp} and equal 
1 in < Cp} with p = pi= C{ph\\ogh\)^. 

3.2 Reduction. Strong and Stronger Magnetic Field 

In this case we have 

(3.34) p > eh-^logh]- 1 , 

(3.35) e = C(p- 1 h\\ogp\)^. 

Then construction of subsection 3.1 works on every its step. Now in a microlocal sense 
l^jl < C{p~ x h\ \ogh\)? (on energy levels below c), and therefore junior terms seem to be 



4 ESTIMATES. 



32 



of magnitude Cp 2 (p~ l h\\ogp\)^ and cannot be skipped; even next junior terms are of 
magnitude Cp?(p~ 1 h\ logp|) 2 which is 0(h) as h\ log/i| 2 < C only. 

In fact, however, as we will see, junior terms will have no impact because a microlocal 
estimate will be replaced by an operator one \rjj\ < C(fi~ 1 h)2 which is sufficient to estimate 
junior terms by Ch as long as p < Coh -1 (strong and superstrong magnetic field). Further, 
for ph > C^h^ 1 (ultrastrong magnetic field) only the lowest Landau level is important. 

4 Estimates. 

4.1 Intermediate Zone. I 

Note Until subsection 4.9 we assume first that magnetic field is intermediate i.e. condition 
(3.1) is fulfilled. 

Let us consider intermediate zone 

(4.1) = {p{ = CmaxOu- 1 , (ph)*) < |£'"| <fh = (M| \ogh\f*}. 
We split it into strips Z p = {p < < 2p} with p in the indicated frames. 

Proposition 4.1. Let condition (3.1) be fulfilled and let Q = Q w be operator with symbol 
supported in Z p , and such that 

(4.2) \d^Q\ < Cp~l Q l 

and ip G C^(B(0,1)). 

Let x e Cg°([-1, -|] U [|, 1]). Then estimate 

(4.3) \F^ h -, TXT (t)TQM < Ch^Tp^ + pT + p,- 1 )^)* 

p z l 

holds as 

(4.4) ~ 2 <T<ep 

p 2 

with arbitrarily large exponent s. 

Proof. Proof repeats those of analysis in subsections 8.1,8.2 of [Ivr6]. First of all, let us 
make 7-partition with respect to (x" , £" , x'") with 

(4.5) 7>C7i- 1 + C r (/r 1 /i|log/i|)3 
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and consider some element ip'{x) of it supported in B(x,j). Let us reduce operator to its 
canonical form in B(x, 27); with no loss of the generality one can assume that this canonical 
form is 

(4.6) /# 2 ^J + A 2 )+ E ti t g»D j D h + V + 0{i) 

l<j<r 2r+l<j<d 

there where fj, gi and V are constant (and we can assume that gi k = Sjf., 0(7) covers 
variations of fj, V as well as unremovable terms and we are discussing only B(x, 27). 
Let us decompose (Tu) into the sum 

(4.7) (Tu)(x,y,t) = £ u a ,p(x",x"',y",y'",t)r a (x')rp(y') 

with T a (x') = r~i \\. v aj (xj/H), h = (/i~ 1 h)'i and Hermite functions v*. 

Then modulo 0(7) in L 2 -norm operator A applied to the terms containing T a (x') in 
this sum becomes 

(4.8) Aa= 7i(2% + l)M+ E ^tf + V 

l<j<r 2r+l<j<d 

which is a temporary notation. 

We strengthen condition (4.5) assuming that 

(4.9) 7 > p 2 + ph + /X" 1 . 

Then we can prove easily that for Q' = QTip'T* the following inequality holds 

(4.10) \F t ^ h - wX T(t)r'Q'i,u a> p\ < C^h'^Tp^x 

h n S / h , s ,p~ l h* 



( / n \s 1 ft u //x ft v 



,p2j>/ v p7 y v 7 2 



^7+ fj(oij ~ otj)\ v 7+ \vhJ2jfj(Pj ~ a j)\ 

with some a = a(r) G Z + r where r'i> = ^ v\ x n=yn )X tr' =y /'i dy" dy'" . 

Here factor ph 1 ~ d p q ^ d comes as a trace norm of Q'u a) p\ the second factor is the sum 
of three terms: (-jf) is due to microhyperbolicity with respect to £ w and rescaling while 



terms (— ) and (^t— ) are operator calculus errors. Finally two last factors in the right- 
hand expression of (4.10) are due to ellipticity of operator hD t — A a . 
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More precisely, let us note that operator A a + 0(7) — r is elliptic in this vicinity as long 
as a 21 = {a G Z +r : | J2j(^ a j + tyfj^h + ^ ~~ T \ < 7} an d we can express there u a> * 
with a G' 21 via u a > ^ with a' G 21. Namely, for a ^ 21 we have 

(4.11) |F t ^ h - lT XT(<)QV/3l < 

vaax\F t ^ h -i T x T (t)Q'u a ',i3\ x (7(1 V(2aj + + ?| + 7) J 

3 

The similar inequality holds for u* t p with /5 ^ 21. Finally, for M aj( g with arbitrary a, [3 2t 

(4.12) |F t _^-i r XT(*)<5X*,/3| < max |F^-i r x T (t)Q'it Q ',/3'| x 

(7(| + 1)/^ + V\+ 7)" 1 )' x (7(1 ^(2/3,. + l)f,ph + V\+ 

j j 

which in combination with microhyperbolicity arguments justifies (4.10). 
Note also that #21 < C^(ph)~ r . Then (4.10) implies 

(4.13) \F t ^ h -i rX T(t)TQ^u\ < C^h-^Tp^x 

h x s , h x S 



((^r+O'+^'H^'" 

CTi - w +'((^)V(A r + ( ^y 



Restoring to original Q, ip (by summation over re-partition) we lose factor r y d . As exponent 
s is large enough the optimal value of 7 honoring (4.9) is 7 = pT + p 2 + because for 
p > p~ l (which we assume) 7 2 > p~ l p 2 T. 

Then we arrive to estimate (4.3). □ 

Proposition 4.1 implies immediately estimate 

(4.14) \F t ^ h -, T {Mt) ~ Xf (t))TQM < CTrV VV + /i^ 1 ) x (^Y 

with x ^ Cq°([—1, 1]) equal 1 on [— |, |], T = e/i" 1 as T satisfies (4.4). 
In turn (4.14) implies 

(4.15) \F t ^ h -, T x T {t)VQM < Ch- i yr 1 P q {^ + P~ X ) x (4)' + C7* 1_ V 



4 ESTIMATES. 



35 



where Cp q h 1 d estimates \F t ^ h -i T Xf (t)TQilm\ due to standard results rescaled. 
Therefore Tauberian arguments immediately imply estimate 

(4.16) TZ 1Q = \T{Q^e) - h- 1 f (F t ^ h -, T x T {t)TQ4>u) dr\ < 

J — oo 

CT-^c/rVVV + z^ 1 ) x (^) S + ^"V) = 

Let q > 2. Integrating the first term in the right-hand expression by — from p\ = 

C(fih)^ to pi we get its value as p = pi; one can see easily that the result will be then less 
than Ch 1 ~ d . This procedure applied to the second term gives 0{h l ~ d ) for sure. 

So for q > 2 in the whole zone {pi < < pi} we get a proper remainder estimate 

(4.17) Tl 2Q = \Y{Q^e) - h~ l f (F t ^rXT n (t)rQ n i;u) dr\ < Ch^ d 



0<n<n 



as Q = J2 < n < n Qn, Qn operators with the symbols supported in {p n < < p n+ i}, 
p n = Tpl pnX pi, and satisfying (4.2), T n = ep n . 
Case q = 1 we analyze later. So, we arrive to 

Proposition 4.2. As q > 2 in frames of proposition 4-1 estimate (4-lT) holds for operator 
Q = Q w with the symbol supported in {pi < < pi} and satisfying {4-2). 

Now let us consider 
(4.18) h~ l [ (F t ^ h -i TX T(t)TQiJu)dT 

J — oo 

where again Q is an operator with the symbol supported in Z p , T G [T ,ep], p G [p*,Pi] 
and I remind x G C£°([— 1, — ~] U [|, 1]). We can rewrite (4.18) as 



(4.19) T-^F^-vXT^rQ^) 



r=0 



with x(t) = it 1 x(t) an d due to proposition 4.1 (4.19) does not exceed 
(4.20) Ch~ d p^ x (A)*. 
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Making a summation with respect to T G [To, ep] we get its value at T = To i.e. Ch~ d p q+2 (ph/ p 2 ) s . 
Then after summation with respect to p G [p\, pi] we get its value at p = p\ i.e. Ch~ d (ph)^ +l . 
Estimating roughly contribution of interior zone {|£ w | < Cp*} by C ' ph l ~ d (p\) q and esti- 
mating contribution of the exterior zone by Ch}~ d due to results of section 2 we arrive to 
estimate 

(4.21) TZ X = \T{tffe) - hT 1 I (F t ^ h -i T x T {t)T^v\ dr\ < Ch^ d + C(ph)i +1 h- d 

J —oo 

with T = T . 

On the other hand, due to the standard results rescaled for q > 2 



(4.22) \hT x j° (F^ h -x T XT{t)ri/>u) dr - Js 



W (x,T)^(x)dx\ < 

C(ph)i +l h~ d + Ch 1 ~ d 



as T = T because with our choice 



Cp p< pi, 

(4 ' 23) £= ^ Cp^(^) s + Chp^\\ogh\ p>fk 

mollification error does not exceed Ch}~ d as well. 
Thus we proved 

Theorem 4.3. Let q > 2 and conditions (0.1) - (0.7) with (I, a) = (2,1), (I, a) = (1,1) 
and (0.24) be fulfilled. Then 

(4.24) n<Ch l - d + Ch- d (ph)i +x ; 
in particular 1Z < Ch l ~ d as 

(4.25) p < p* q = Ch~&. 

Note that p* 2 = Ch~i, pi = Ch~i etc. 
Now we can assume that 

(4.26) p* q < p < eh~ l \\ogh\' 1 . 
In this case we will use the following corollary of (4.16): 
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Proposition 4.4. Let q > 2 and (4.26) hold and Q = Q(hD"'), Q = Q(g") be operator 
with symbol supported in Z p , pi < p < p x , ij) £ C^°(B(0, 1)). Then 

(4.27) K lQ = \T(Qe) - h' 1 f (F t ^ h - lT XT(t)rQ^u) dr\ < Ch^p^ 1 

J —oo 

as T = ep and 

(4.28) TZ 0Q = \F(Qe) - h' 1 I (F t ^ h -x T x To (t)rQi>u) dr\ < 

J — oo 

ch l - d p^ l + ch- d (^_ypi. 

pi i 

with T = ep. 

Note that the right-hand expression of (4.27) integrated by y does not exceed Ch 1 ~ d . 
4.2 Intermediate Zone. II 

Now we can apply the technique of the previous subsection to calculate expression 

(4.29) f (F t ^ h -, T <p T [t)TQ^u) dr 

J —oo 

with <p = x- Let us consider the difference 

(4.30) f (F t ^ h -i T <j) T (t)TQil){u - u'j) dr 

J —oo 

between (4.29) and the same expression for perturbed operator A' = A + 0(p~ 1 ). Then 

(4.30) is equal to Iq+ I\ + ■ • • + /«, where Iq is defined by (4.30) with 4> = x an d T replaced by 
T = ep" 1 and /„ are defined by (4.30) with 4> = x an d T = 2 n T , n — 1, . . . , n, n = [T/T ] . 

Applying arguments of the previous subsection and noting that an operator norm of 
1 tA _ & ih 1 tA'^j does not exceed Cp~ 1 Th~ 1 as \t\ x T, we conclude that /„ does not 
exceed (4.20) multiplied by Cp~ 1 Th~ 1 as n > 1, i.e. 

(4.31) | f (f^-^xt^TQ^u - u')) dr\ < Ch~ d p^ 2 x (A;)' x p~ l hr l T. 

J —oo P ■*- 
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After summation with respect to n we get the same expression with / = To: 

(4.32) | j° (F t ^ h -i T (x T (t) - XT (t))rg^( M - u')) dr\ < C^hT^ p q+2 {^) S '. 

Note that the right-hand expression in (4.32), integrated by — with p from p* to 1, does 

not exceed its value as p = pi i.e. Cp~ 2 h~ d ~ 1 (ph) 3 ^ = Ch 1 ~ d (ph) !L 2~ . 

Therefore for q > 2 and Q supported in intermediate zone expression (4.30) is equal to 
Iq modulo 0{h}~ d ). However, for I we can apply standard theory rescaled and replace it 
by the difference of the Weyl expressions for A and A' (we will do it in section 5). So far 
we proved 

Proposition 4.5. In frames of proposition 4-4 

(4.33) |r(Q#) - h- 1 ! (F t ^ h - wXT {t)TQiJu - 

J — oo 

*Ufc-irX? (*)W(« - «°)) dr\ < Ch^ d 
where u° is a Schwartz kernel of e lh 1<A ° , A = T*A°T, 

(4.34) A od ^ p 2 J2(^ aJkW (k+h 2 D 3 g^D k + aZ 

j,k<r j,k>2r+l 

is the main part of A. 



4.3 Interior zone. I 

We will return to the intermediate zone (for q = 
zone 



1) later; now let us consider an interior 



(4.35) fi = {|ri<Pi = WH 

still assuming that q > 2 and (4.26) holds. 

As q > 2 one can see easily that the perturbation of the magnitude 0(/i _1 ) there leads 
to an error Oip^p^h^) = 0(h}- d ) in £% w . 

Thus as q > 2 instead of A we can consider in Qq U Qi reduced operator .4° with coef- 
ficients a jkw = a^ix'^x'^p^hD"), gi kw = g^ k {x\x'\p~ 1 hD"), a% = a Q (x", x'\ p^hD"), 
(J = p^hDj + ixj and a jfc = unless 3m : j,k e m; this will lead to an error 0{h 1 ~ d ) in 
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IZ2 q. However we provided different arguments for interior zone fio and intermediate zone 
Qi to justify this replacement. 

If fj had constant multiplicities we would have even 

(4.36) A° = p 2 J2f 3 {h 2 D 2 + p 2 x 2 -ph)+h 2 E D]g^Dl + aZ 

j<r j,k>2r+l 

and decomposing u into (4.7) series we would arrive to the family of operators 

(4.37) A a = h 2 D^ kw D k + W:, W a = a + Y,2afihf J 

j,k>2r+l j<r 

and we would be able to apply without any significant modifications [Brlvr], [Ivr5]. 

However, we need to consider more general operator (4.34). Consider operator Q = Q w 
with symbol supported in Z p with p < p\. 

Repeating arguments of the proof of proposition 4.1 we arrive to estimate (4.13) where 
Q' = TQtp'T*, ip' is an element of 7-admissible partition in (x",x"',£"), 7 > ph. 

Restoring to Q we arrive to estimate 

(4.38) \F t ^ TXT {t)YQM < Ch-«Tp«j((JLy + (A)' + 

_ 1 

As q > 2 and p > p* we can take 7 = ph and p > po > p 9 ; remaining inner core 
^00 = {W"\ < Po} "will be considered separately. Then both h/(p r y) = l/(pp) and p~ 1 h/'j = 
l/(p 3 h) do not exceed h 5 and we arrive to estimate 

(4.39) \F t ^ h - wXT {t)TQM < Cph^Tp^^-Y 
which in turn implies 

(4.40) \F^ h -i T x T (t)FQM < Cph l ~ d p q x A = C/i^-y- 2 , 

P 2 

Therefore, contribution of Z p to the remainder estimate is 0(ph 2 ~ d p q ~ 3 ); integrated over 
f it gives 0(ph 2 ~ d ) = 0(h l - d ) as q > 4, 0(ph 2 ~ d \ \ogh\) = 0{h}- d ) as q = 3, 0(ph 2 - d p^) 
as q = 2. 

On the other hand, contribution of the inner core f2 o = < Po} to the remainder 

estimate is O(ph 1 ~ d p q ). 
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Therefore we will get remainder estimate O^h 1 dj rfih 1 d p q ^) as q 
hpQ 1 )) as q = 2 and we need to pick-up 

(4.41) p = max(/i3 ; fjT?} 

finally arriving to the remainder estimate 0(h l ~ d ) as q > 3 and 0(/i 1 - d + ///iI- <i ) as g = z 
n 

Thus we proved 

Proposition 4.6. Let q>2, (4.26) be fulfilled and Q = Q(hD"'), Q = Q(£"') be operator 
with symbol supported in {|£ w | < Cp\, ip G Cq°(B(0, 1)). 
Then 

'0 q > 3, 



(4.42) n 2Q <ch 1 - d + 



Ch3~ d p q = 2, 



where x e C£°([-l, 1]), x = 1 on [-|, |], <3 = 5-) <»<nQ"> are operators with symbols 
supported in {\p n < W'\ < 2p n } as n = 1, . . . ,N, and in < 2p } os j = 0, p n = 2 n po, 
n = [logpi/poJ +1, T n = ep n (n > I), T = ep _1 . 

This statement remains true with u replaced by u° which is the Schwartz kernel of 



gih 1 tAo 



4.4 q = 1 in Generic settings. Intermediate zone 

Let us repeat arguments of subsection 4.1 first as q — 1, 

(4.43) /Ts| log/i|~3 < p < e/T 1 ! \ogh\~ 1 , 

(4.44) x p > p* = C(p/i)5 + Cp" 1 . 

The first problem is that the factor (p 2 + p~ x ) in the right-hand expression of (4.14) 
is too large as p < Ch~*; furthermore even if there are no unremovable cubic terms, one 
should take factor (p 2 + p _1 /i| log h\) instead because we need to take Cp" 1 /^ log h\ scale 
with respect to x<i. I handled this problem in [Ivr6] as d = 3 and the idea from there (may 
be combined with the idea used in subsection 2.3) is applicable now. 

The factor problem is rather simple. 

Let us assume first that g^ h , Fj^ are constant (and in particular there no unremovable cu- 
bic terms). Then this problem arrises only as p 3 < h\ \ogh\ which implies p < h~^\ log/i|i. 

11 Similarly we would get 0(h 1 ~ d + fihi~ d ) as q = 1 but there is a lot of other things to consider. 
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Let us consider a 7-covering with 7 = p s and consider separately two types of elements: 

(i) elements with |VV| > £ = |log/i|~ 2 . If |c?a. d V| > ( we have that £<j-shift for time 
T e [T , p 51 ] is of magnitude no less than (T and then since logarithmic uncertainty principle 
is fulfilled 

fj,-^ x e xp" 2 C > C7i|log/i| 

we derive that xt^Q^u is negligible; otherwise |Vj_y| > ( and similar arguments work for 
(x", f ")-shifts. 

(ii) elements with |VV| < (. We can take on them 7 = C(p~ l h) a | log/i| which is larger 
than Cf) -1 /i| log/i| and replace the term 7 by C(C7 + 7^ \ogh\~ a ) which does not exceed 
C\xh. In this case one can take T\ = epQ~ l . 

Similar analysis combined with ideas of subsection 2.3 works in the general case as long 
as there are no unremovable cubic terms. 

If there are unremovable terms of this type one can replace decomposition (4.7) by a 
similar decomposition with T a orthonormal eigenfunctions of operator 

(4.45) a = A + /U" 1 Yl PjkmLjLkLm + CoP^Al 

j,k,m<2r 

with 

(4.46) Aq= f At' 1 1' 1 ■ •'•})• 

L2J-1 = fihDj, L 2 j = Xj and Tp are replaced by their complex conjugates). We added the 
last term in (4.45) to make operator non-negative and self-adjoint; this term is 0(fi~ 2 ) on 
u and thus it is included in the approximation error estimate 0(p~ l \ log/i|~°") anyway. We 
got operator (4.46) by rescaling x' \—> Dj 1— > pDj. 

For this operator a in L 2 (W) the standard semiclassical eigenvalue asymptotics holds 

(4.47) n(r, h) = Weyl + 0{h l ~ r ), Weyl x tT r as K -> +0 
and in particular 

(4.48) n(r + h,h)-n{r,h)=0{h 1 ~ r ) 

where 'Weyl" means Weyl expression for a and we plug h = ph. 

Then the same modifications hold for estimates (4. 14), (4. 15) as well and then Tauberian 
arguments improve estimate (4.16) to 

(4.49) Ki Q <c(rWx (4r + C7^)|log P r 2 
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as p < h 5 1 because T\ > e\p\ logp| 2 now. 

Integrating the right-hand expression over ^ from p\ to pi we get 0(/i 1_d ) then. However 
for h 6 ' 1 < p < h" 1 ] \ogh\~ 1 we take only T\ = e^p and recover estimate 0{h 1 ~ d \ \ogh\). 

On the other hand, under microhyperbolicity condition for p < eh~ l \ \ogh\~ 1 we can 
take Ti = e\ and recover estimate 0(h l ~ d ) again. 

Thus we have proven 

Proposition 4.7. Let q = l, (I, a) = (1,2). Then 

(i) Estimate (4-24) holds. In particular for p < Ch~z sharp remainder estimate 1Z < 
Ch l ~ d holds. 

(ii) in frames of proposition 4-2 TI-2Q defined by (4 -IT) does not exceed h l ~ d \ \ogh\ where 
again T n = ep n . 

(Hi) If either p < h 5 " 1 or fj have constant multiplicities or microhyperbolicity condition 
is fulfilled at level 0, 1Z 2 q defined by (^.1?) does not exceed Ch l ~ d where T n = p n \ logp„| 2 . 

(iv) Furthermore, if Q is supported in {|£ w | > p} and equal 1 in {|£ w | > p} then IZoq 
defined by (4.28) does not exceed Ch l ~ d + Ch~ d p(ph/ p 2 ) s . 



4.5 q = 1 in Generic settings. Interior zone 

In contrast to the case q > 2 skipping 0(p~ l ) terms is not now generally justified. Because of 
this we do not try now arguments of subsection 4.2 and we will use arguments of subsection 
4.3 without skipping such terms first, assuming that 

(4.50) ChT^ < p < eh-^logh^ 1 , 

(4.51) CpT 1 < \£ d \~P<P*i = W i 

Starting from (4.38) with 7 = ph due to the above arguments we arrive to 

(4.52) \F t ^ h -i T (x T (t) ~ Xn(t))TQilm\ < 

Cph^p(T (4-) s + T(A)- + T(^y 
V V To PI 7 

as T = 4 < T < Ti = ep and thus 

(4.53) \F t ^ T x T {t)TQM < Cph l - d p(- 2 + T(— ) s + T(^M . 

\p PI 7^ / 
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Therefore due to Tauberian arguments 
(4.54) \r(Q^e) - h' 1 f (F t ^ h - w x T {t)TQ^u) dr\ < 

J— oo 

CyT-W-tpf^ + T( -Y + T(^) s ' 
\p z PJ 7 



cp/^p(/^ 3 +(^r+(^)' 

V jot t 



*P7' T 

After integration the right-hand expression over ^ with p ranging from p , 

(4.55) p = Ch~s > Cp" 1 as p > p*, 

to p]; we get 

-y) +cph^ d P {^ 

-fTf' ) p=po Y p=pi 



(4.56) Cph l ~ d p(hp- 3 + (-Y) + Cp/^V^V 



where all terms excluding Cph 2 d p 2 = Cpht d do not exceed Ch 1 d . 

Also, contribution of the inner core {|^| < Po} does not exceed Cph 1 ~ d p = Cph^~ d as 
well. 

So, the the final estimate is Ch 1 ~ d + Cplv$~ d . Thus we arrive to 

Proposition 4.8. Let q = 1, condition (4-50) be fulfilled and Q = Q(hD"'), Q = Q(i"') be 
operator with symbol supported in < Cp~i, G C^°(B(0, 1)). 
Then 

(4.57) K 2Q < Ch 1 " 11 + Cphl~ d 

where IZ2Q is defined by (^.17) with Q = ^2 0<n< nQ n {hD'") , Q n are operators with symbols 
supported in {\p n < W'\ < 2p n } as n = 1, . . . , n, and in < 2p } as n = 0, p n = 2 n pQ, 
n = [logpi/poJ + 1, T n = ep n (n > I), T = ep _1 . 
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4.6 q = 1,2 under microhyperbolicity condition 

We repeat arguments of subsection 4.3, but we will use microhyperbolicity at some moment. 
Namely, (4.10) implies that 

(4.58) \F t ^ h -, T xr(t)TQ'4m\ < C^h'- d Tp'x 

7 s / I tr(A(x',x"',Z") + t~f)~ s \\ dx"dx"'d^" 

JB{z,2-y) 

where 5(2,7) ^ suppQ', -4 is an operator reduced and considered as an operator in 
L 2 (W,C) and ||.|| is an operator norm there. Then the same estimate (with an integral 
over -8(0, 1)) holds as Q' is replaced by Q and due to microhyperbolicity the last integral 
is of magnitude ^{ph)~ r as 7 > p 2 ; even condition 7 > p _1 is not needed anymore if we 
consider operator (4.44) from the beginning 12 . 
So, 

(4.59) \F MXT {t)TQM < Ch^Tp* x (( A)' + (A)' + (AA)*) 

as h/p 2 < T < ep and picking up 7 = p 2 , p > po = h^~ s > p" 1 we rewrite the right- 
hand expression as Ch~ d Tp q+2 (— — Y . Note that we changed here definition of p an d the 

v p l T' 

corresponding notion of the inner core. 
Then (4.59) implies that 

(4.60) \F t _ h - dT x T {t)TQM < Ch^pi 

as long as T < ep, p > p . 

However, due to the microhyperbolicity condition F t _ h ~d T XT{t)TQipu is negligible as 
long Ch\ \ogh\e~ 1 < T < e and T = ep satisfies this inequality as long as 

(4.61) p>Cp , £> C/ip _1 |log/i|. 

So, (4.61) yields that (4.60) holds for T < e. 

12 As q > 2 we can skip all the 0(// _1 ) terms. As q = 1 we preserve all the terms which do not contain 
factor hDd, we skip all such terms containing (hDd) 2 because this leads to the approximation error p 2 ^ 1 
and then we remove as in [Ivr6] all such terms containing hDd exactly in power 1. 
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Then we can apply Tauberian arguments giving estimate with right-hand expression 
Ch 1 ~ d p q . This last estimate hold for p < pi as well. Integration of this expression over y 
gives remainder estimate 0(h l ~ d ). 

Consider now contribution of the inner core {|£ w | < po}- In this zone we get estimate 

(4.62) \F t ^ h - dT x T (t)TQM < Ch~ d p q + 3 , 

as T = p and p ==po there; then microhyperbolicity arguments push it to T = e under 
condition (4.61) and then Tauberian arguments lead us to contribution to the remainder 
not exceeding Ch~ d pl +3 which is 0(/i 1_d ) as q > 1, 5 small enough. 

Thus under microhyperbolicity condition remainder estimate is 0(h 1 ~ d ) even as q = 1,2 
but we need to look at an approximation error; one can prove easily that for smallest 
e == Cp _1 h\ log/i| satisfying (4.61) it is also 0(h l ~ d ). Thus we have 

Proposition 4.9. Under microhyperbolicity condition TZ±q does not exceed Ch 1 ~ d where 
Hiq is given by (4- Iff) with T = e. 

4.7 q = l,2 and constant multiplicities of fj 

Assume that fj have constant multiplicities. Then in both interior and intermediate zones 
after reduction and skipping 0(/i _1 ) terms which are due to the third resonances and all 
smaller terms 13 instead of "matrix" operator we get a family of "scalar" operators 

(4.63) A a (x", x'", pT x hD'\ hD'") = 

J2 (hD 3 )g jk (x'\ x'", pL- l hD'")(hD k ) + V a (x", x'\ pT^D'") 

j,k<r 

and we can study them separately. Let us apply standard rescaling procedure: for each 
index a we introduce functions 7 = 7 a and g = g a hy 

(4.64) 7 = infjr : 0(r) > p 2 + \V a (x" ,x"' ^'%$(r)r~ l > |W a (xV",£")l} + 7, 

(4.65) = 0( 7 )3, 

(4.66) £ = #(7)3, Q^ = Ch =^ q = ft,*+2 1 log pi\~ 20+27 , 7 = /iiT2 1 \ogh\~& 

where we assume that (1, 1) ^ (/, a) ^ (2, 0). Note that 7 2 > CpT x h\ log 

13 The case of q = 1 where there are such terms will be considered in details in the next subsection 
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(i) Consider first elements with 7 > C7. Calculating contribution of each element 
of (7; p)-partition with respect to (x", x'", £"') we can apply elliptic arguments unless 
€ Q < p < Cq\ in the latter case we can apply microhyperbolic arguments; then contribution 
of each element to the remainder estimate does not exceed C p r h r ~ d+1 g q ~ 1/ ~f d ~ l in the latter 
case and 

(p + 0)7 

in the former one. 

Then the contribution of each such element W of 7-partition with respect to (x", x'", £') 14 
to the remainder also does not exceed C ' p r h r ~ d+1 g q ~ 1/ -f d ~ 1 . 

Then the total contribution of such elements to the remainder estimate for given a does 
not exceed 

c J ffhr-^g^-^dz 

with z = (x",x"',£") and after summation over all indices we get 

(4.67) C ^' hr ~ d+1 j et^dz. 

(ii) On the other hand, contribution of each element of 7-partition with 7 < C7 does 
not exceed C p r h T ~ d ^ d Q q and their total contribution does not exceed (4.67) as well. 

Since for given (z, p) 

(4.68) #{a G Z +r : ep 2 <el< cp 2 } < C{ph)- r (p 2 + fih) 
expression (4.67) does not exceed 

(4.69) Cfih 2 ~ d J erSr 1 dz + chl ~ d J p'lipy 1 d p 

where 71 = min a 7 Q and Qi corresponds to it, j(p) is defined from equation -#(7) = p. 

One can see easily that the second term does not exceed Ch}~ d , while the first term 
does not exceed 

Cph 2 / 7 2 |log7| 2 c?7 



which in turn does not exceed Cph l ~ d g q unless q = 2, / = 2. In the latter case o < and 
the first term does not exceed Cfih 2 ~ d \ logh\~ a . In both cases we arrive to the remainder 
estimate 

(4.70) TZ 2Q < Ch 1 ^ + Cph^+^l \ogn\~Tm 

14 More precisely, of W x < Cp\}. 
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which is exactly remainder estimate (0.40). I remind that it is Cph 2 ~ d as (I, a) = (2,0). 
Here T = j/g. 

Thus we have proven 

Proposition 4.10. Let q = 1,2, p > p*, fj have constant multiplicities and A contain 
no unremovable cubic terms 15 . Let Q = Q(hD"'), Q = Q(£ w ) be operator with symbol 
supported in < Cpi}, i> G C^(B(0, 1)). 

Then estimate {4-70) holds with IZ2Q given by (4-17) with Q = J2o<n<n Qn{hD"'), Q n 
are operators with symbols supported in {^7 n < 7 < 27„} as n = 1, . . . , n, and in {7 < 27} 

as n = 0, N = Llog7j + 1, 7n = 2"", T n = C T ^| log 7n |-f fa > 1/, T = e^" 1 . 

On the other hand, if q = 1 skipping unremovable terms in ^4 comes with a price 
of 0{^h 1 ~ d + fi~2hz~ d ) error. If we would prefer not to remove them, we cannot treat 
separate equations. However we can do as well as before assuming that 

(4.71) fj = const, j = l,...,r 

which by no means excludes variable g^ k , or unremovable cubic terms. 
We can introduce then 

(4.72) 7 = infjr : 0(r) > /o 2 ,0(r)r -1 > \VV(x", x"\ + 7, 

and q by (4.65) and 7 by (4.66). Then we can apply our standard arguments as long as 

(4.73) p 2 + I W| > Cp~ x 
which will be fulfilled automatically for 7 > C7 as long as 

_ 1 _ 21 2a 

(4.74) p > C/i 2 p>Ch w\\ogh\w 

with the right-hand expression exceeding h~i for sure. So, under condition (4.74) estimate 
(4.70) holds. On the other hand, if condition (4.74) is violated, we need also consider 
elements with violated (4.73), in particular, |VV| < Cp" 1 and 

(4.75) 7 < 7 = Cp~T\ log/i|T , 

15 May be we just dropped them in the interior zone causing an error 0(h 1 ~ d ) if either q = 2 or q = 1 
and each of these terms contained factor pT^hDd and an error 0(^h 1 ~ d ) in the generic case as q = 1. 
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Instead of 7 partition we consider max(7, ^-partition and on elements with 7 = 7 let us 
introduce 

(4.76) 7 ' = 7/i p 2 + 7 / , g ' = p + p\ 

(4.77) 7' = ^373/13 = ft 3// 31 |10g/l|M, p = fJt, 3-y 3/2,3 

Then we arrive to the remainder estimate Ch 1 ~ d + C ph 2 ~ d ^'~ l : 

Proposition 4.11. Let q — 1, p > /^ and condition (^.71) 6e fulfilled. Then 

(i) Under condition (4- 74) estimate (4-70) holds. 

(ii) If condition (4-74) is violated then estimate 

(4.78) TZ 2Q < Ch l ~ d + C/i/i 1 "^^! log/xpT^ + C/i^/if- d | log/x|-a 
holds. 

Remark 4.12. (i) Condition (4.74) means exactly that the third term in the right-hand 
expression of (4.78) does not exceed the second one. 

(ii) If condition (4.71) is violated one needs to take in account that as the main part is 
of the form Q + Ei<j< r fji^) + h 2 D 2 ) + ^ ^ fi = £"), V = V(x", x d , £"), 
its gradient with respect to x",Xd,C' depends on the localization with respect to {(p 2 x 2 + 
h 2 £ 2 ),j — 1, . . . ,r} and there is at least one obstacle to it: due to logarithmic uncertainty 
principle these quantities are defined with Cph\\ogh\ precision. Therefore one must take 
7 /_1 | log7| _<T | > Cp,h\ logh\ in the arguments leading to proposition 4.11 which deteriorates 
estimate (4.70). It does not look worth needed efforts. 

(iii) Another approach would be to drop unremovable cubic terms as p > q paying 
the price of 0{h 2 ~ d Q~ l ) for this and apply the same arguments as above. Then remainder 
estimate is Cph^~ d g3'j3 where 7 = £?T|log/i|T and one should optimize the answer with 
respect to q which is not pT^ anymore. 

4.8 Improved remainder estimates 

There are peculiar cases when without microhyperbolicity one can improve generic remain- 
der estimate TZi < Ch 1 ^ + Cph 1+ i~ d or estimates (4.70), (4.78). It happens only for 
r > 2. 

Note that the factor estimating the number of indices one needs to take in account is 
actually (n(r, ph) — n(r — Cp 2 , ///i)) where n(r, ph) is an eigenvalue counting function for 
operator a in L 2 (W, C) calculated at any point (x", x"',£") at suppQ while so far we used 
estimate C(p 2 + p,h)(p:h)~ r for this difference. 
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Therefore one can replace term Cph 1 d g q we have in frames of proposition 4.10 by 



Pi 



(4.79) Cp r h r ~ d+1 J j (n(x,r,fj.h) - n(x, r - Cp 2 , ph)) p^ 1 ^ 1 dp dx+ 

Cp r h r ~ d J (n(x, r, ph) — n(x, r — Cg 2 ,ph))g Q dx 

with 7 = gT| log p | T , g = /i i+2 | log/i| ; + 2 where we already returned to coordinates x. 

Also, in the generic settings one can replace term Cph 1+ 3~ d by expression (4.79) with 
7 = p 2 and g — hs. 

Consider two examples 

Example 4.13. Consider eigenvalue counting function no = no(r, h) of operator ao = Ao in 
L 2 (M. r , C); since Ao depends on (x", x", £"), so does n and we can always return to variable 
x by map \Ev Since Ao contains only quadratic terms one can calculate n explicitly 

(4.80) n (r, h) = #{a G Z +r : ^(2o^ + l)f)h + V < r). 

j 

If all fi , . . . , f r are comeasurable then 

n (r, K) - n (r -X,h)< Cti~ r (\ + h) 
is the best possible estimate (A > 0); otherwise from 

(4.81) n (r, K) - n (r - A, h) = ^ #{«' G Z +r " 1 : £ ^(2a, + l)H G 

n>0 j>2 ■' 1 

[ /r i (r _ V ) _ 2 („ + 1) - /^A, f;\r -V)- 2{n + 1))} 

one can derive a better estimate Ch~ r {\ + ^(fr)) with u(h) = Ch K , k > 1 depending on 
Diophantine properties of {/2//1, • • • , fr/ fi} (see [Har]). 



Example 4.14. Let us assume that 



ICi 



(4.82) rank{V^,...,V^} >k-1 Vz. 
Ji Ji 

16 Actually wc need a quantitative version of this condition and only at points where W is a linear 
combination of V/i, . . . , V/ r - 
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It is easy to recover an integrated version of (4.81), namely (0.33) with u(h) = Ch K . Re- 
ally, one needs just to consider (^-admissible partition in (ah, x)-space with scaling function 

C = c|V(r-V-E i (2a i + l)/ i n)|. 

Furthermore, these arguments show that for n the left-hand expression of (0.33) does 
not exceed Ch l ~ r (\ + h K + 

Then we arrive to 

Proposition 4.15. Let condition (0.33) hold for n which is an eigenvalue counting func- 
tion for Ao- Let fi* < /i < eh' 1 ] \ogh\~ 1 . 

(i) Let either q = 2 or q = 1 and A contain no unremovable cubic terms 17 . Then 

(4.83) K 2Q < Ch 1 ^ + Cu{fih)hi- d ; 

(ii) Assuming instead that n which is an eigenvalue counting function for a, satisfies 
(0.33), we need neither assumption "A contains no unremovable cubic terms" in (i) nor 
(0.33) for n ; 

(Hi) Let q = 1 and A contain unremovable cubic terms and fi > h~z \ log | " 2 . Then 

(4.84) K 2 ,q < Ch 1 ^ + Cu(fih)h^ d + Cfi^h^. 

Proof. Statements (i) and (ii) follow directly from above arguments. In (iii) one needs to 
notice that (0.33) for n implies 

(4.85) J (n(x, r, h) - n(x, r - A, h)) dx < C(X + ^ + u(K))h- r . 

□ 

Remark 4.16. (i) If q — 1 and (0.33) holds then dropping unremovable cubic terms leads 
to an error not exceeding 

(4.86) Cu(fih)fi~h~ d + Cfi~h~ d . 

(ii) As fi > h~3 last terms in the right-hand expression of (4.84) and in (4.86) do not 
exceed Ch 1 ^' 1 . 

Similarly one can prove 



1 Unremovable cubic terms containing hD c i are allowed. 
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Proposition 4.17. Let q = 1,2 and f) have constant multiplicities. Let condition (0.33) 
hold for n . Let fi* < fi < eh~ 1 \ log 

(%) Lei either q = 2 or q = 1 and A contain no unremovable cubic terms. Then 

(4.87) ft 2 Q < C/i 1 ^ + • | log W ; 

(m^ Let q = 1, condition (4-71) be fulfilled, A contain unremovable cubic terms, estimate 
(0.33) hold for n as well and fi > h~z | log h\~% . Then estimate (4-84) holds as (4-74) ^ s 
fulfilled and estimate 

(4.88) TZ 2Q < Ch 1 '* + Cv(^h)(h- d+ i^\logii\-m +Cfi~'^rh-l- d \logh\-^ 
holds as (4-74) fails; 

(Hi) Let q = 1, condition (4-71) be fulfilled, A contain unremovable cubic terms and 
/i > h~z 1 \ogh\~2 . Then estimate 

(4.89) TZ 2Q < Ch x ~ d + C(v(/ih) + ■ h~ d+ ^\ \ogfi\~^) . 
holds as (4-74) is fulfilled and estimate 

(4.90) n 2Q <ch 1 ~ d + 

C(u(fih) + fi' 1 ) ■ [h~ d+ TT2\ log/i|"^) + Cfi-^rh-^~ d \ \ogh\-^ 

holds as (4-74) fails. 

4.9 Strong Magnetic Field 
4.9.1 Generic settings 

Consider now strong magnetic field case eh~ 1 \ \ogh\~ 1 < fi < eh' 1 . I remind that reduction 

to the canonical form was done with e = C(/i _1 /?,| log/i|) 2 rather than CyU -1 ; furthermore, 

exterior zone is empty now and f>\ = C rather than f>\ = C(/i/i| log/i|) 2 . 

As q > 2 throwing away 0(p~ 1 ) terms brings an 0(h 1 ~~ d ) error of £ MW . Then all the 
arguments of subsections 4.1-4.3 remain true and remainder estimate will be 0(h 1 ~ d ) as 
q > 3, 0(jj,h^~ d ) as q = 2 and 0(fih~3~ d ) as q — 1; in the latter case contribution of the 
intermediate zone is estimated by Ch 1 ~ d \ log/i|. So, in the generic settings previous results 
hold without modihc&tions: 

Proposition 4.18. Propositions 4-1, 4-%> 4-4~4-8 and theorem 4-3 remain valid for strong 
magnetic field as well. 

Thus in what follows we need to consider cases q = 1, 2 only. 
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4.9.2 Microhyperbolic case 

We must distinguish two cases: 

(i) Microhyperbolicity direction £ does not depend on r. In this case all the arguments 
of subsections 4.1-4.4 remain true and the remainder estimate is 0(/i 1_d ). 

In particular, the case of fj having constant multiplicities is covered. Really, in this case 
microhyperbolicity condition means exactly that at each point x 

(4.91) £ r,V{h = 0, V, TJ > =► V, Tj = 0; 

l<j<ri 

then there exists vector £ = £(x) such that 

(4.92) (4v(i))>0 Vj 

which exactly means microhyperbolicity in direction I. 

(ii) Microhyperbolicity condition is fulfilled with £ depending on r. Then as in [Ivr7] 
one should take 

(4.93) e = C7i|log% c ^ |log/l1 

and while remainder estimate will be still 0(h l ~ d ), approximation error does not exceed 

Ch~ d (^ l h\ log/i|)5 e c ^ |logh| | \ogh\-° 
which is 0{h l ~ d ) as / > 1 and jih < e' with small enough e'; as I = 1 we arrive to 

(4.94) Hi < C/ i 1 - d e c ^ |log ' l| | log/il 1 -*. 
So, we have proven 

Proposition 4.19. In i/ie strong magnetic field case under microhyperbolicity condition 
TZi does not exceed Ch l ~ d where TZi is given by (4-16) with T — e, Q — I provided either 
microhyperbolicity direction does not depend on r or I > 1. Otherwise estimate (4-94) 
holds. 

4.9.3 q = 1,2 and fj have constant multiplicities 

Assuming that there are no cubic terms we get a family of separate operators A a and all 
the previous arguments hold bringing us remainder estimate (4.70). On the other hand, as 
q — 1 dropping unremovable cubic terms produces an error C^h l ~ d which is less than the 
right-hand expression of (4.70) unless I = 2, a £ (—4,0]. In this case assuming (4.71) one 
can prove easily (4.70) as well. So we have 
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Proposition 4.20. (i) In the strong magnetic field as fj have constant multiplicities es- 
timate (4-70) holds unless q — 1, I — 2, a G (—4,0] and there are unremovable cubic 
terms. 

(ii) In this exceptional case should include extra term C^h l ~ d in the right-hand expres- 
sion. 

4.9.4 Improvement without microhyperbolicity 

Furthermore, under condition (0.33) for no we get remainder estimate (4.83) in frames of 
proposition 4.15(i) which is a very subtle improvement now. 

Moreover, if fj have constant multiplicities we can recover (4.87) in frames of proposition 
4.17 (i). 

Proposition 4.21. In the strong magnetic field propositions 4-15, 4- 17 remain valid. 

4.10 Superstrong Magnetic Field 

Consider now the superstrong magnetic field case eh^ 1 < fi < Ch" 1 . Then the same 
arguments as before but in simplified form hold and one can prove easily 

Proposition 4.22. For superstrong magnetic field for 

(i) In the generic settings remainder estimate IZ2 < Ch x ~ d holds as q > 3, remainder 
estimate 7Z 2 < Cfi0~ d = 0(h3~ d ) holds as q = 2 and remainder estimate 7I2 < C[ihz~ d = 
0(hz~ d ) holds as q — 1. 

(ii) If q = 1,2 and microhyperbolicity condition 0.2 with t independent on r is fulfilled 
then remainder estimate 1Z\ < Ch 1 ~ d holds. 

(Hi) If q = 1,2 and fj have constant multiplicities then remainder estimate 
TZ 2 < C/T d+ A I log /i I ~2$+2j holds. 

Note that we can ignore cubic terms (which are 0(h) now). 

4.11 Ultrastrong Magnetic Field 

The last and the easiest case to consider is fi > Ch" 1 . Then we arrive to a single operator 
(4.95) A = h 2 D j g jk {x",x"',fi~ 1 hD")D k + a Q {x",x"',fi- 1 hD") 

2r+l<j,k<d 

where I remind that 

a = jyo^ , W = V-^2f jt ih 

j 
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belongs to C l,cr uniformly. 

Surely we get a system, but all other components u a $ (with (a, (3) ^ 0) could be 
expressed via u ,o- 

Then the principal part of asymptotics (as \W\ < c) is of magnitude fi r h r ~ d . For the 
remainder estimate we have 

Proposition 4.23. For ultrastrong magnetic field 

(i) As q > 3 remainder estimate IZ2 < Cfi r h r+1 ~ d holds. 

(ii) As q = 1, 2 and microhyperbolicity condition 

(4.96) \W\ + \VW\>e 

is fulfilled then remainder estimate IZi < C // h r+l ~~ d holds. 
(Hi) As q = 1,2 remainder estimate 



Cn r h r+ ^- d \hgh\~& q = 2, 

Cn r h r+ TT2- d \\ gh\-TT2 q = i 



(4.97) TZ 2 < 



holds. 



On the other hand, mollification and approximation error in operator is 0(h) and if 
either q > 2 or microhyperbolicity condition holds mollification and approximation error in 
Weyl expression is 0(h) multiplied by fi r h r ~ d i.e. fi r h r+1 ~ d . 

On the other hand if q — 1 and there is no microhyperbolicity condition mollification 
and approximation error is O(h^) multiplied by n r h r ~ d i.e. fi r h r+ 2~ d . 

In both cases mollification and approximation error does not exceed the remainder 
estimate. 

5 Asymptotics 

Now we need to derive more explicit asymptotics rather than formula (2.27) or its clones: 

h 1 / F t ^ h -i T x T (t)r(Q x i/; x u) dr, 
J —00 

(5.1) h- 1 / F t ^ h -i T XT(t)n^ y Ql)dr 

J —00 

with different T where ip y = i()(y), Q y = Q(y,hD y ) and Q l means the dual operator; we 
write operators acting with respect to y on the right from a function. Both expressions 
coincide due to the trace property. 
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5.1 Weak magnetic field redone 

Weak magnetic field case is when ji < h 5 ~ l and we derived asymptotics with the principal 
part given (5.1) with Q = I , T = and with some remainder estimate derived without 
canonical form reduction. 

Note that the mollified operator A is microhyperbolic in direction (£, d^) due to condition 
(0.7). Further, in zone {|£ w | > e'} this operator is a differential operator since we take there 
e = Ch\ \ogh\ which does not depend on |£ w |. 

On the other hand, in zone {|£ w | < |e} operator A is microhyperbolic in direction 
+ (£"?<%') an d it is a differential operator with respect to x',x" since mollification 
parameter e = Cp~ l h\ \ogh\ does not depend on 

Due to the standard results rescaled then F t _ h -i T XTif)^ iuip y ) is negligible as |r| < e' 
and T G [To, To], T = Ch\ log 7i| , T = e/i" 1 . Therefore, due to (2.26) we can take T = T 
in (5.1) with only 0(h s ) difference. It happens in the proof of theorems 0.3, 0.5 and 4.3. 

Then we can launch successive approximation method with the unperturbed operator 

(5.2) A = Y,Pf9' k Pk + V, 

3,k 

9 jk = 9 jk (y), V = V{y), h = hD k - f jV k (y)-^{d k V){y){x k -y k ) 

j 

and plug it into (5.1). We follow [Brlvr]. In what we get then the first term does not exceed 
CTh~ l ~ d = Ch~ d \ \ogh\. Perturbation can be written as 

with /i-differential operators and since each factor (xj—yj) according to transformations 
of [Brlvr] leads to extra factor T in the estimate, one can see easily that each next term in 
the successive approximations gains factor Ch~ 1 T(T 2 + /iT 3 ) < Ch\ logh\ 18 . 

Therefore only first two terms should be considered. The first term will give exactly 
(5.1) with u(x,y,t) replaced by u(x,y,t) and it will be exactly J £ MW (y,r)ip(y) dy with 
0(h s ) error 19 . 

The second term consists of two parts; one of them is generated by the perturbation 

(5.3) Ri=J2 P A d m9 jk )(y)^m - Vm)Pk + ^2(d m V)(y)(x m - y m ) 

j,k,m m 

18 I remind that factor Th^ 1 comes from Duhamel principle. 

19 After we replaced u by u we can replace T by oo with 0(h s ) error since then (2.26) holds for any 
T > ep with 0{h s T- s ) error and it will be exactly / £ MW (y, r)\jj(y) dy. 
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and it is obviously 0, while the second part is generated by 

(5.4) R2 = J2 fa - p j)9 3k Pk + Pfg ]k {P k - P k ) + (Pj - i) W h iJ\ - h)) 

and is obviously 0(h~ d ~ 2 fiT 3 ) = O^h 1 ^). 

Therefore in this case we can replace (2.27) by J S MW (y, r)ip(y) dy without deterioration 
of the remainder estimate 20 

► This takes care of Theorems 0.3 and 4.3. Also under microhyperbolicity condition 0.1 it 
takes care of Theorem 0.5. 



5.2 Intermediate magnetic field. Decomposition 

Now let u consider case of intermediate magnetic field pf < \i < eh~ 1 \ logh\~ l . Then 
there is a standard zone where we use already derived estimate for 1Ziq with T = Tq; this 
standard zone contains exterior zone > pi} but could be wider. Now in this standard 
zone we can replace T by T = Ch\ \ogh\. 

The unexpectedly difficult (as q — 1) problem is to join asymptotics derived in the 
standard zone with cut-off operator / — TQT* and in the remaining non-standard zone 
where we employ cut-off operator TQT* . To tackle it better let us rewrite the formula we 
derived for the answer: 

(5-5) h~ l Y, / ^WvXT^rHvQ^dT 

where Q n are appropriate elements of the partition of / and T n are already chosen. The 
asymptotics with this principal part and different remainder estimates were derived in 
section 4. Actually one can replace here T n by larger values without affecting remainder 
estimate, but since we have not proven this we will do a bit differently. Considering some 
term here with Q = Q n and T = T n one can rewrite it as 

(5-6) h" 1 ( Yl <P m ,L m {r)[F t ^ h - lT Xf{t)T{u^ y Ql))dr 

J -co < m<oo 

where supp0 o C [-2,2], L = Cpi, supp0 m C [-2, -1], L m = 2 m L as m > 1 and 4> m ,L m 
form partition of 1: J2o< m< oo 4>m,L m (r) = 1 as r < 0. 

20 We need also to replace £ MW by £ MW but estimate of the approximation error is already done; in what 
follows we will not distinguish between £ MW and £ MW keeping in mind that we always must include an 
approximation error estimate in our final statement. 
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Note that one can rewrite each term with m > in (5.6) as 



(5.7) m ,U^A)(^(*)r(«^)) 



<P = X, T = T. 

t=o 



Let us consider (5.7) with <p = x an d T > hL^ . Using arguments of subsection 4.1 can 
estimate such term by 

(5.8) ^-VLiTx(A)^ 

1 i-'ra 

where factors L m , T are due to integrals over r, t respectively and another factor L m is due 
to calculation of the number of contributing indices a; we assume that 

(5.9) L m >C(^ l + p 2 ). 

The same Ch~ d p q L m estimate holds for (5.7) terms with tp = x and TL m = 1 as well. Then 
the sum with respect to t-partition results in Ch~ d p q L m , and then the sum with respect to 
r partition results in Cp q h~ d because maxL m x 1 here; terms with L m > C are less than 
Ch s L^ due to standard ellipticity arguments. 

Let us employ method of successive approximations described in subsection 5.1 and plug 
it into (5.7). The estimate Cp q h~ d for the first term we have as well. 

A perturbation with factor (xj — y,j) will result in either factor CT 2 h~ l if commuting 
with A or A or in in factor CTp~ l if commuting with Q or, what is equivalent, in factors 
ftX~ 2 and hL^p~ l respectively thus resulting in 

(s.io) ch^L- 1 x ctf-y-* x (-^-Y 

respectively. 

Note that the commutator with Q n (and thus the second of these expressions) appears 
in the sum only as T > eT; as for T < eT commutators from adjacent elements compensate 
one another; I remind that T = T n . 

Therefore taking the sum with respect to t-partition and then with respect to m we get 



(5.11) Ch>- d p q L»\ Ch l - d p q - 1 x (_ y 

respectively. We can always take L = pi since T n > /i" 1 for sure and get O^h 1 ^). As 
q > 2 we can do even better but it does not matter. Then both terms in (5.11) do not 
exceed Ch}~ d and remain this way after we sum over partition in zone {|£ w | < pi}- 
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So, in terms with m > we can completely ignore partition and consider only the first 
term of successive approximations, Further, after this we can take T arbitrarily large (thus 
we can take T = +oo); then after easy calculations we get 

(5.12) f [(l-^ Lo (r))d T S MW (x,r)ij(x)dx. 

J —oo J 

So we proved estimate 



(5-13) \h~ x f {l-(f>o, Lo (r))J2( F ^rXT n 

J —oo „ 



(t)(r«V w g* ) )dr 



i - <p 0yLo (T))d T £ MW {x, t) ij}{x) dx\ < Ch 1 -*. 
Now we need to consider terms with m = i.e. 

(5.14) h- 1 f o , Lo (r) (^-vXrWrHQy) rfr ' T = T - 

J —oo 

Instead we consider the same expression but with T = Tq = efi^ 1 and consider a correction 
later. We remember that T n > To in the intermediate and interior zones. So we are 
considering (5.14) with T = T . 

Now we want to replace T by a lesser value. We could replace it by T = Ch\ \ogh\ 
with a negligible error before but it is not the case anymore since as a result of mollification 
and transformation symbol of operator satisfies the only 

(5.15) \d«a\ < C^ M \]ogp 1 \- ff + C 

in the intermediate and interior zones. However then we can take there 

(5.16) T = Cp^h\ log h\ <f . 

Now let us apply the same successive approximations to calculate (5.14) with T = T . 

Thus perturbation with the factor (xj — yj) results in T 2 h~ l = Pi 2 h\ \ogh\ 2 factor and 
extra pj~ 2 |log/i| 2 is well absorbed by p q+1 at least as q > 2, p < log h\~ 2 . One can 
consider q = 1 and the exceptional case either using rescaling arguments thus punishing 
T > h/pi by factor {h/{p\T)) or just by taking the two-terms approximation instead of 
the one-term and proving that the second term is identically 0. 

Furthermore, perturbation containing factor p(xj — Uj)(xk — Uk) will get an extra factor 
pT z h~ l = ph 2 \ logh\ 3 pi 3 which should be treated in the same way. I leave details to the 
reader. 
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So, we can again consider only the first term of successive approximations but in contrast 
to the previous part we cannot tend T = T n to infinity in this approximation term unless 
in the exterior zone. Then the answer will be similar to given in (5.12) 

(5.17) h- 1 I <p , Lo (r)(F t _ th ^ T XT (t)r(u4j y ))dr = 

J — oo 

f j (0(r) + Lu(r))d T £ MW (x, r) ^{x) dx 

J — oo J 

with 

(5.18) w{r) = Toh- 1 J° <P(T>)(T'n( {T ' ~ h T)T ° ) dr' - 0(r), 
where <p = o ,l o - 

Here co(r) ^ because we cannot replace T by +oo. Instead we replaced back T by T 
since operator A has the same propagation with respect to x properties as A and therefore 
even with T y instead of T the difference is negligible. 

Let us define 



(5.19) 



^ w (x, o) = ^r^ 2 ^- v £ / fi---fr^L Q(n dr, 



(5.20) n a (x,T) = {£:V + J2( 2a i + 1 )f* h + £ + ^ < t}, 

i 2r+l<j,fe<d 

where gi k with j, = 2r + 1, . . . , d denote coefficients after we reduce operator A to its 
canonical form 

(5.21) £ /,(/r/;j • fx)) • £ g^D^ + V, 

l<j<r 2r+l<j,k<d 

g'" = <l<u (//'%:. 2r+] „. 

Then £ MW = + Plugging this into (5. 12), (5. 17) and adding we get in our 

calculation candidate to the final answer 

(5.22) J S MW {x,0)^{x)dx + J° J u{T)d T S MW {x,T)^{x)dx = 

£ MW (x,0)ij(x)dx+ f [ u(r)d T £™ w (x,T)ij(x)dx 
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modulo negligible error where Q is supported in < 2pi} and equal 1 in {|£ w < 2pi}. 
Now it is a time to remind about correction 



because we replaced T n by T in the intermediate and interior zones. Thus the following 
statement is proven: 

Proposition 5.1. Let conditions of proposition 4-4 be fulfilled. Then expression (5.5) with 
summation over all zones modulo (^(/i 1 "^) is equal to 



where U = T*ur*\ ip = T*tpT, Q n = T*Q n T and here we can take Q n = Q n (hD"') 
elements of partition in zone < 2pi}, Q = Q n . 

Now what we need is to estimate the sum of the second and the third terms in the 
right-hand expression of (5.24) (but not separately). 

5.3 Decomposition. II 

Let us consider remaining inexplicit terms in the the the sum in expression (5.24) 



where we renamed 4>o,l into 0(r) (so, in (5.18) <p is no more 0o,l o ) an d again for the sake 
of simplicity we denote T n by T and skip index n. 

Let us apply the same decomposition technique to 4>{r) as we did to 1: 0(r) = 
Tl,m>o 4>m,L m ij) with lesser L m than before and rewrite (5.25) as the sum of (5.7)-type 
terms 



(5.23) 




(5.24) 




(5.25) 




(5.26) 
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Replacing (xf(t) ~ Xf {t)) by Xt(^) we find that each obtained term with m > does 
not exceed (5.8) again; however, there is a difference: after summation with respect to T 
ranging from To now to T we get that expression (5.26) does not exceed 

(5.27) Cpih- d L m x A s 

J-'m 

as long L > ph and m > 0. After summation with respect to L m it gives us Cpp q h x ~ A . 
It would enable us to set L = ph if not condition L > p 2 + + ph we need. 
Let us assume that yU > h~^\ then we can take 

(5.28) L = Cp 2 + Cfih; 

analysis in the case q = I, h~3 < p < h~z will be done later with or without microhyper- 
bolicity condition. 

Let us employ the same method of successive approximations with unperturbed oper- 
ator A = Ao(y", y'", pT x hD" , hD'") considered as operator with operator valued symbol in 
L 2 (M. d ). Then the perturbation is 

(5.29) - vjM^ t/', y'"' »~ lhD "> hD "') + h ( x "> f"> t/ H , ^ lhD "> hD "') 

r+l<j<d 

with b = Oi/j- 1 ). 

Consider the second term of approximation plugged into (5.26) with m > and trans- 
form it in our usual way. Factor (xj — yj) in the first part of perturbation (5.26) should 
commute with either Q ± bringing thus factor pT 2 h~ x as 2r + 1 < j < d or smaller factor 
p~ x T 2 h~ x as r + 1 < j < 2r or with Q bringing factor p~ l T . However again commuting 
with Q in the final run (after summation with respect to n) should be considered only as 
T > T = T n and since T > ehp~ 2 due to section 4 the first factor is larger anyway. 

Thus we get from the first part of perturbation terms estimated by (5.8) multiplied by 
one these factors: so we get 

(5.30) ! Ch- d p"Lx (^-Y x pT 2 h'\ 
(5.30) 2 ChT d p q L x x ^T 2 h~\ 
(5.30)3 Ch-'pOLx^YxTp- 1 
respectively because in our analysis T > e/i" 1 , L > ph. 
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On the other hand the second part of the perturbation (5.29) brings factor C p, x Th 1 
or equivalently CpT x L~ x thus giving us 

(5.31) C IT 1 (fifC* x (A)*. 

Summation of (5.30)i-(5.30)3, (5.31) with respect to L, T results in the same expressions 
but with the minimal possible values, i.e. T = pT x and L = p 2 + fih; so the result for each 
of the expressions (5.30)i_ 3 does not exceed 

(5.32) C^p^h- d x(^-) s , 

p l + fih 

while the result for (5.31) does not exceed 

(5.33, c„-yv- x (-J^y. 

Integration over -£ brings instead of (5. 32), (5. 33) their values as p = p\\ the first one is 

always 0(h l ~ d ) while the second one is 0(/i 1_d ) as q > 2 only and is 0(fi~^h^~ d ) as g = 1. 
We can do better than this but marginally. We will do it later. 
So, replacing U by its one-term approximation 

(5.34) U = -ih~ x G+b{x - y)8(t) + %hr x Q-b[x - y)8(t) 

in each term (5.26) with m > brings the total error 0{h 1 ~ d ) unless q — 1 and A contains 
(unremovable) cubic terms; in the latter case the total error is 

So, let us plug U into (5.26). I remind that the symbol of i/j is given by (3.28). But then 
the remainder estimate 0(/i -1 | \ogh\~ 1 ) there multiplied by p\p q h~ d estimate (I remind that 
extra factor p\ comes from support of <fi) results in 0(h 1 ~ d ). 

Further, terms containing factors r]J with j = 1, . . . , 2r result in because these terms 
applied to v (x", x"')T a (x') result in p:\a- p\=i v /3( x "y x"')T a (x') while Aq "honors" \a\. 
Finally, terms containing r]J with j — 2r + 1, . . . , d also result in since we can take Q n 
even with respect to and Aq is even as well. 

So we need to plug only ip instead of ip. We can also replace T n by any larger value, 
say infinity. But then using change of coordinates \l/o we g e t exactly the second term of 
(5.24) with the opposite sign and with uj defined by (5.18) for <p = <fi mt L m - Then we arrive 
to (5.24) with L = ph + p 2 (which now depends on n). 

We also need to look at the case m — 0. So far we considered both intermediate and 
interior zones. Consider now intermediate zone only. Then all above arguments remain 
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true with the only difference that factor (-J?f) S is now due to estimate (4.3) rather than 
T-mollification. Therefore 

(5.35) \J2 f I u;(r)d T 8^(x,r)^(x)dx+ 

h-'Yl [ ^o.io( T ) ( F t^h-w{xT n (t) - XT (t))^U'ijjyQ t ny S ) dr\ 



with summation over n such that Q n belongs to intermediate zone does not exceed Ch l ~ d 
if either q > 2 or q = 1 and A contains no (unremovable) cubic terms and Cp~^h^~ d if 
q = 1 and A contains (unremovable) cubic terms. 
Then from proposition 5.1 we arrive to 

Proposition 5.2. Let conditions of proposition 4-4 be fulfilled and leth~^ < p < h" 1 ] \ogh\~ 1 . 

Then one can rewrite expression (5.5) as (5.24) with L = Cph and <f>, uj defined 
accordingly and one should take Q n = Q n {hD"') elements of partition in interior zone 
< Pi}; Q — J2 n o^nd the error is 

1. Oih 1 ^) if either q > 2 or q = 1 and A contains no (unremovable) cubic terms; 

2. 0(p~^h^~ d ) if q = 1 and A contains (unremovable) cubic terms. 

I remind that the only cubic terms which matter are those which do not contain hDj 
with j > 2r + 1. 



5.4 Interior zone. Generic settings 

Now in view of proposition 5.2 we need to consider only interior zone < pi}- In 

the generic setting there is no point to take L < ph because we cannot improve estimate 
0(ph l ~ d ) for the number of indices making contributions as r runs an interval of length 
L . Because of this (as m = 0) one should replace (5.30)i_ 3 by 

(5.36)! Cpp^h 1 ^ x min((^) s , l) x P T 2 h~\ 

(5.36) 2 Cpp^h 1 ^ x min((^) s , l) x pT x T 2 h~ x , 

(5.36) 3 Cp,p q h l - d x min((^) s , l) x Tp~ l 
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respectively where factor (y^) s as T > h/p 2 is due to estimate (4.3). Furthermore, (5.36)3 
one should count with T > eT n x p only. 

After summation with respect to T we get instead of (5.36)i_2 their values as T = h/p 2 
and instead of (5.36)3 its value as T = p i.e. 



(5.37)i_ a Cuffl-*}?-*, Cp"- 4 h 2 ~ d , C^h^i^Y 

p A 



respectively. After summation with respect to p > h$ we get 0{h l d ) as q > 3, 0(p,hs d ) 
as q = 2 and 0{p,hi~ d ) as g = 1. 

On the other hand instead of (5.31) we get 



(5.38) Ciip q h l ~ d x min((^) s , l) x pi^Th' 1 

and summation with respect to T results in C p q ~ 2 h l ~ d . 

Summation of this expression with respect to p should be taken over p > max(/i"2, hs^j 

only; it results in C^/i 1 " 6 *) as q > 3, 0(/i 1 ~ d | \ogp, 2 h\) as q = 2 and 0(p^h 1 ~ d ) as g = 1. 

To get rid of factor | log /i 2 /i| as q = 2 one should consider two-term approximation, then 
(5.38) will be replaced by 



(5.38)* Cpip q b>- d x min((^) s , l) x pT 2 T 2 h 



-2 



which results in C pT x p~~ 2 h x ~ d after summation with respect to T and in 0{h l ~ d ) after 
summation with respect to p. 

On the other hand, consider the second part of perturbation (5.29); one can rewrite it as 
b(y" ,y" f , fi~ 1 hD" , hD'"). Really, since this part contains unremovable cubic terms, symbol 
b belongs to JF 1 ' 1 and thus freezing x" = y", x'" = y'" one makes an error which could be 
accommodated in the first part of (5.29) and treated correspondingly. 

However, the second term in the approximation of U is 

(5.39) U' = -ihQ + BQ + 5(x - y)8(t) + ihQ-BQ-5{x - y)5(t). 

and plugging ip instead of z/> results in trace equal because operators Aq and thus Q ± 
transform v(x",x"')T a (x') into the sum of vp(x" ,x"')T p(x') with \fl — a\ — 0, 2 while B 
transforms v(x", x"')T a (x') into the sum of vp(x", x"')Tp(x') with \/3 — a\ = 1,3 if we 
consider terms bjkmVjVkVm with of j < 2r, k < 2r, m < 2r. 

Terms with j < 2r, k < 2r, m > 2r + 1 etc and j > 2r + 1, k > 2r + 1, m > 2r + 1 are 
odd with respect to £ w while everything else is even and thus disappear after integration 
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with respect to Furthermore, plugging (tp — i/jq) instead of ip produces extra factor p" 1 
in the estimate which is more than enough to compensate logarithm. 
Now we need to estimate contribution of zone {|£ w | < Po} with 

1 2 

• po = h~i if either there are no unremovable cubic terms or p > h~s and 

1 2 

• Po — if there are o unremovable cubic terms and p < h~z. 

However it follows from section 4 that we can consider this zone as a single partition 
element and set T n = To = ep~* here thus bringing extra term C , poT ~ 1 /i 1 ~ rf = Cp\ph x ~ d 
to the remainder estimate; the latter expression is equal to 0{h l ~ d ) as q > 3, 0(ph3~ d ) 
as g = 2; as g = 1 this expression is equal to Cphs~ d if either p > or there are 
no unremovable cubic terms in A and it is equal to Cp^h l ~ d if p < /i~3 and there are 
unremovable cubic terms in A. 

So with this error one can replace U by U in (5.25). Further, we can replace ij; by ipo 
there with 0(/i 1_d ) error. 

Now we can preserve remainder increasing T n to oo. Really, consider 

(5.40) h- 1 f (f>oM( F t^rXT(t)rU^OyQl)dT 

J — oo 

with T > Chp~ 2 ; it does not exceed 

(5.41) Ch^xphxp^x (A)' 

p z i 

exactly as would be for U but now we can take any T > Chp~ 2 because A is a constant 
coefficient operator-valued operator. Summation with respect to T > T n results in value of 

(5.41) with T = T n which is p in these settings i.e. we get 

(5.42) Ch^^ + phjpox^y 

and after summation with respect to p we get its value as p = which is C phi +l ~ d . 

But then after obvious calculations the second and the new third term in formula (5.24) 
just cancel one another and we are left with just the first term. 

So we had proven 

Proposition 5.3. Let conditions of proposition 4-4 be fulfilled and leth~^ < p < /i" 1 ] log/i|~ 
Then one can rewrite expression (5.5) as J £ MW (x, 0) ip(x) dx with an error 
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1. 0(h l ~ d ) asq> 3; 

2. 0(h l ~ d + phl~ d ) asq = 2; 

3. 0{pht~ d ) as q = 1 and either p, > h~% or A contains no (unremovable) cubic terms; 

as q = 1, fi > h 3 q = 1 and A contains (unremovable) cubic terms. 

Remark 5.4. (i) When we say unremovable cubic terms we mean terms bjkm^j^lk^lm with 
j < 2r, k < 2r, m < 2r only. Really, any other term bears extra factor p and if A contains 
only these terms one can take po = Chi and the above analysis of subsections 5.3, 5.4 
produces error 0(/i 1_d ) (we get rid off factor | \ogp 2 h\ in the same way as we did it for 
q = 2 in the proof of proposition 5.3. 

(ii) In subsection 5.5 we show that under microhyperbolicity conditions unremovable 
cubic terms do not cause any corrections even as q — 1. 

(iii) In subsection 5.6 we derive asymptotics with a better remainder estimate. 

(iv) In subsection 5.7 we derive asymptotics with a correction due to unremovable cubic 
terms. 

► Proposition 5.3 together with remainder estimate of proposition 4.4 and with mollification 
error estimate imply theorem 0.8. 

5.5 Calculations under microhyperbolicity condition 

If either q = 2 or q = 1 and there are no unremovable cubic terms, everything is easy: 
under microhyperbolicity condition one can replace in the estimates (5.36)i_3 factor \ih by 
factor p 2 as p > p = CTis; then instead of (5.37) i_ 3 we will get 

(5.43)!_ 3 cvrt 1 ^, CKy- 2 /^, Cp q+2 h~ d {^-) s 

v p 

respectively which will result after integration over dp/ p in 0(h l ~ d ) in each of them 21 . 

Further, contribution of the inner core is estimated in section 4 by CT^ 1 pg +2 h 1 ~ d = 
0(ph 2 ~ d ) as we take T n = Tq there. So we can again replace U by U and then ip by ipo 
and after this T n by oo. Then after obvious calculations we estimate expression (5.26) by 
Ch 1 ~ d and thus (5.24) will be reduced to its first term modulo 0(h l ~ d ). 

21 Sure in the first expression as q = 1 we get C/i 1_d | \ogh\ but we consider two-term approximation as 
above and also we replace ip by ipo; then the error term results in 0(/i 1-d ) while the second term (5.39) 
resulting in 0. 
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On the other hand, factor ph is replaced by p 2 in (5.38) as well thus resulting after 
summation with respect to T in Cp~ A p q h~ d ; after summation with respect to p we get 
Cp~^h^~ d as q — 1. Even use of two-term approximation gives us (5.38)* multiplied by 
p 2 /(ph) and after summation with respect to T and then by p we get Cp~ 2 h~s~ d which is 
0(h 1 ~ d ) only for p > h~~s. 

However, if A contains no unremovable cubic terms in the sense of remark 5.4(i), one 

should replace there factor p" x T by p~ x pT\ thus we get Cp~ 1 Tp 4 h~ d min^ (7^2) S , lj in 

both cases instead of (5.38). Then after summation with respect to T we get Cp~ 1 p 2 h~ d 
and after integration over dp/ p we get C p~ x p\ 2 h~ d = Ch 1 ~ d . 

Also, we need to cover also the case h~z | \ogh\~z < p < h~* (lesser p are covered 
by weak magnetic field case due to microhyperbolicity) but one can weaken assumption 
p > h~z if A contains no unremovable cubic terms. 

So, we can plug U instead of U (and ipQ instead of ijj) and thus we have proven 

Proposition 5.5. In frames of proposition 4-4 f or I 1 — h~*\ \ogh\~^ under microhyperbol- 
icity condition 0.1 one can rewrite (5.5) as J £ MW (x, 0)ip(x) dx with 0(h l ~ d ) error provided 
either q = 2 or q = 1 and A contains no unremovable cubic terms (in the sense of remark 
5.4(i)). 

However, we want to get rid off this extra assumption and thus we want to consider case 
u q = 1, A contains unremovable cubic terms". To do this we must reexamine intermediate 
zone as well. 

Our main idea is to estimate decay of F t _>h-i T XT{t)FijjoQU more accurately, using mi- 
crohyperbolicity conditions it was done in subsections 4.1-4.3, [Ivr7]. 

Proposition 5.6. In frames of proposition 4-4 f or t 1 — h~*\ \ogh\~^ under microhyperbol- 
icity condition 0.1 

(5.44) \F t ^ TX T{t)Ti>QU\ < CT p^h^^) 1 ^ log(A) |- x min((A)* ; i) 
as po = h~^ < p < pi, |t| < e. 

Proof. Repeats ones of propositions 4.1, 4.2 of [Ivr7] and is based on ^-approximation of A 
with hT- 1 <?]< ChT~ l \\ogh\. □ 

After we improved our estimate by extra factor (^)' 1 \ log(^) \~ a we can do more precise 
estimates. 
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Obviously, 

(5.45) /T 1 f <f> 0>Lo (T)(F t ^ h -i TXT (t)'ipQnU)dT = 

with x(t) = t^xif) and (f>' derivative of O ; I remind that (f>o(0) = 1. Here L Q could be 
larger than p 2 to accommodate p < p~2. Then proposition 5.6 implies that for q = 1 the 
right-hand expression of (5.45) does not exceed 

(5-46) CA- d (^) / " 1 |log(^)|-xmin((A)^i). 

Integrating (5.46) over dT/T from T — Tq to T — oo we get this expression as T = To 

1 112 

and integrating over dp/ p from p = hz to p = p\ as /i~2 1 log 7i| _ 2 < p < h~s we get this 
expression as p = Pq i.e. 

(5.47) C(p/i)^|log/ir CT /^ d 
which does not exceed approximation error 

(5.48) Cp~ l p 1 1 log h\- a h- d = Cp^- l hh~ d \ log/i|5- CT 

as p < Zi~2. So, in this case we can plug U instead of U (and ipo instead of 

Thus it remains to consider case /i~2 < p, < ft, - 3. Consider (5.46) with T = p _1 and 
the last factor equal to 1. Integrating it with respect to y from p = h~s to p = p~2, we 
get its value as p = p~^ which is 



(5.49) Cp-5(p/ l ) i - 1 |log/i|- CT / i - d 

and which does not exceed approximation error as I < § and C/i 1 "^ as Z > §. These 
arguments estimate contribution of the inner core as well. Therefore we estimated properly 
contribution of inner zone < Cp~2}. 

Consider now p ranging from p~2 to p\. We also apply the above arguments as T > 
fih > resulting after integration over dT/T in 

(5.50) Cp*p}~ l \\ogh\-°{\) s h- d . 

pp A 
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As To < T < ph we apply two-term approximation method with the second term 
resulting in and with an error estimate 



(5.51) 




which can be easily proven by arguments of propositions 4.1,4.2, 5.5 of [Ivr7]. After inte- 
gration over dT/T we get 

(5.52) Cp 2i -V 2 | logp|- CT /T d x min((^) s , l). 

Finally, integrating (5.50) over dp/ p in the indicated limits results in its value at the lower 
limit p = p~z; this result does not exceed p~^~ l \ \ogh\~ a h~ d which again does not ex- 
ceed approximation error (5.48). This is also true for (5.52) as I < |; as I — | we get 
Cp~ 2 \ \ogh\~ a h x | log p 2 h\ which also does not exceed approximation error. 

So again we can replace U by U and ip by ipQ. Further, after this we can replace T n = e 
by T n = oo with a negligible error. Thus we have proven 

Proposition 5.7. Let q — 1. Then in frames of proposition 4-4 f or f 1 — h~^\logh\~^ 
under microhyperbolicity condition 0.1 one can rewrite (5.5) as J S MW (x, 0)ip(x) dx with 
0{h}~ d + pk~ l hk~ d \ log/i|5~ CT ) error. 

► Proposition 5.7 together with remainder estimate of proposition 4.9 and with mollification 
error estimate imply theorem 0.6. 

5.6 q = 1, 2 and constant multiplicities of fj 

If we drop all junior terms we get a family of scalar operators and everything is easy. 
However covering case q = 2 this would not cover case q = 1 even if unremovable cubic 
terms contain hDd- So we need more delicate arguments. 

Note that instead of Q n (hD"') we have now Q n (x", x"\ pr x hD")Q/ m (hD'") where Q' m are 
elements of p-partition, Q n are elements of 7-partition with 7, q introduced by (4.64)-(4.65). 

Consider first elements with q > Cg and p > eg. 

Then one should multiply (4.36) by 7 d x g 2 /j where factor 7^ comes from Q n and factor 
g 2 /'-/ (which is greater than p' 1 ) comes from estimate of | V x ",x'",£'"A\. 
Then (5.37)i_3 become 

(5.53) !_3 Cpg^h 2 -^- 1 , Cg^h 2 -^' 1 , C pg^h 1 ^ m.m{{^)\ l) 7 d - 1 



5 ASYMPTOTICS 



70 



respectively and the second and the third expressions do not exceed the first one. 

On the other hand, as p < eg, g > Cg we can rescale problem to microhyperbolic with 
x i— > 27, h 1 — > h/(g / ~f), p~ l h 1— > p~ l h'-f~ 2 . Then the total contribution of elements with fixed 
n and different m does not exceed (5.53)i again. So, (5.53)i estimates the contribution of 
element Q n . 

After summation over partition (5.53) 1 gives 

(5.54) ph 2 ~ d J g q - l 1 - 1 dx"dx"'di" 

which does not exceed the right hand expression in (0.41) (or (4.70)). 

Finally the same estimate remains true as g x g because now we estimate the contribu- 
tion to each of the second and the third term of (5.5) rather than to their sum. 

Similarly, if unremovable cubic terms contain factor hDd then we get that the contribu- 
tion of Q n instead of (5.48) is equal to 

(5.55) C^-^Vmin^" 1 ^ 1 , l) < Ch l - d g q ~ l 1 d 

and the total contribution of all elements does not exceed Ch l ~ d because then p scales as 
HQ 22 ■ 

On the other hand, if there are unremovable cubic terms without factor hDd then (5.55) 
is replaced by 

(5.56) C y u/i 1 -y 7 d min(/i- 1 £- 2 , ^ < c^i-y-y 

and summation results in Cpih 1 ~ d and Ch l ~ d as q — 1, 2 respectively. 
Thus we arrive to 

Proposition 5.8. Let q = 1,2 and fj have constant multiplicities. 

(i) Let either q = 2 or q = 1 and there are no unremovable cubic terms. Then with error 
equal to the right-hand expression of (4-69) one can rewrite (5.5) as J S MW (x, 0)i/j(x) dx; 

(ii) Let q = 1 and there are unremovable cubic terms. Then with error equal to the 
right-hand expression of (4-70) plus Cp^h l ~ d one can rewrite (5.5) as J S MW (x, O)ijj(x) dx. 

► Proposition 5.8 together with remainder estimate of proposition 4.10 and with mollifica- 
tion error estimate imply theorem 0.11. 



22 



Parameter /1 1 h scales independently from [i and h because we deal with the reduced form of operator. 
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5.7 Improved asymptotics without microhyperbolicity 

Now it is easy to prove 

Proposition 5.9. In frames of proposition 4-15(i) one can rewrite (5.5) as 

J S MW (x, O)ijj(x) dx with an error not exceeding the right-hand expressions of (0.34), (0.36) 

for q = 2, 1 respectively. 

► Proposition 5.9 together with remainder estimate of proposition 4.15 and with mollifica- 
tion error estimate imply theorem 0.10(i),(iii). 

A bit more sophistication requires 

Proposition 5.10. In frames of proposition 4-17(i) one can rewrite (5.5) as 

j ^mw^ O)ip(x) dx with an error not exceeding the right-hand expressions of (0.41), (0-42) 

as q = 2,1 respectively; 

Proof. One needs just to replace ip by ipo, then consider two term approximations (the 
second term will result in in the end), and then replace T n by 00. In the intermediate 
zone estimates repeat those of the generic case. 

In the interior zone we need to consider integral (5.54) and notice that contribution in 
it of all elements with c~ 1 t < g < ct, t < (fih)? does not exceed 

Ch}- d Q q - 1 1 - 1 (Q 2 + V(^h)) 

with q defined by given 7 and thus this integral does not exceed 

Ch l ~ d J Q q+1 1 - l dx + Cv(^ih)h l - d J gf-^dx 

and plugging g — 72 1 log7| - ^ and replacing dx by 7~ 1 c?7 we find that the first term does 
not exceed Ch l ~ d provided either q = 2 or q — 1, (I, a) y (1, 1) while the second integral 
does not exceed 

CVGu/i^-V" 1 ^" 1 x Cv(fih)e q h- d 
unless q = 2, / = 2 in which case the second term does not exceed Cu(fih)h 1 ~ d \ \ogh\~ cr . □ 

► Proposition 5.10 together with remainder estimate of proposition 4.17 and with mollifi- 
cation error estimate imply theorem 0.12. 
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5.8 q = 1: asymptotics with correction 

We want to improve remainder estimates in asymptotics as q = 1 and A has non-removable 
cubic terms. To do this we need to include some correction in £ MW . 

First of all, if both \ an d 4> are proper functions 23 and L > Cph\ log/i|, we can rewrite 
(5.5) in (5.24) form with 0(h s ) error because (5.7) with <px = Xt and with m > is 
negligible as TL m >Ch \ \ogh\. 

To calculate the last term in (5.24) we apply the same successive approximation method 
as before but with unperturbed operator A = Ai(y", y'", p~ x hD" , hD'") where in A\ we 
include unremovable cubic terms 

(5.57) M = fi 2 (bjkmVjVkVmY 

j,k,m,<2r 

but do not include similar sum with at least one of the indices j, k, m equal d . 

Then the last term in perturbation (5.29) will contain at least one factor rjd and thus 
we get extra factor p in the corresponding estimates. Then instead of unwanted term 
Cpph}~ d x (pp 2 )~ 2 calculated as p = max(/i~2,p ) in an estimate we get Cpph 1 ' 11 x (pp)~ 2 
calculated as p = max(/i"^, p ), which results in 0{h}~ d )] here we take the same two- 
term approximation. We still need to consider the second term to avoid extra log h factor. 
However, we need to consider it only as b contains exactly one factor hDd because otherwise 
we gain one more extra factor p. 

We also can replace ip by tpo with impunity because the we get an error not exceeding 
CLopp~ l h~ d = Cp\ log h\h l ~ d which does not exceed Ch 1 ~ d as p < h~ 1 \ \ogh\~ z and is 
much less than the remainder estimate otherwise (actually we can replace in the above 
estimate Lq by ph but we do not need it). But then the second term in this approximation 
will be odd with respect to ^ and result in in the final answer. 

So, all estimates of propositions 5.3, 5.8, 5.9 related to the case "q = 1 and A contains 
no unremovable cubic terms" hold for an error arising if we replace in the last term of 
(5.24) U by U which is the first term in our new approximations and also replace ip by ipQ. 
Therefore in comparison with our previous calculations we get an extra term 

(5.58) h^J2 f 0o,L o (r)(^ h - lr (xTjt)-Xf„W)r((?7-Wo J/ Q^))rfr 

n J -°° 

where U is an old first term in the approximation. 
23 In [Brlvr] sense. 

24 We also include C^ 2 A^ to have operator semibounded from below. 
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This is a correction term in question but we want to rewrite it in more explicit form. First 
of all we can replace here T n by +00 (using the same arguments as before) and also replace 
f o by f = Ch\logh\ because (F^ h -i T XT(£)r(£A/> yQ^)) is negligible as T e [T,T ], 
|t| < CL for both U and U. 

However if we consider 

(5.59) h^Y, f fo,L (T)(F^ h - iT Mt)r{(U-UWo v Q t ny ))dT 

and apply method of successive approximation to calculate it (considering Ao as an unper- 
turbed operator and A\ as a perturbed one, we can get easily an error estimate Cp\ log h\ 2 h l ~ 
for it which is again less than the remainder estimate we are looking for. 

Thus modulo term not exceeding remainder estimate we can rewrite (5.58) as the same 
expression but with xt replaced by 1 and with Xf replaced by i.e. 

(5.60) h^J2 f <f>o, Lo (T)(F t ^ h -i T r(U -U^oyQ^dr 
which in turn is exactly equal to J £^?(x, 0)ip(x) dx with 

(5.61) = (2tt)-^^-V /(I " 0o,L o (r)) x 




x, t - z 2 ) - n (x, r - z 2 )) Q(z) dz x f x ■ ■ ■ f ry /g. 



Thus we arrive to 

Proposition 5.11. Let q = 1, h~^\ log/i| - ^ < h~ 1 \\ogh\~ 1 and A contain unremovable 
cubic terms. LetTZ 3 be the difference between expression (5.5) and 

J(£ M ^(x,0)+S^(x,0))^x)dx 

where E^^ is defined by (5.61) with Lq = Coph\ log h\ with large enough constant Co and 
a proper symbol Q equal 1 in {\z\ < Cpi} and no and n are eigenvalue counting functions 
for operators ao and a respectively. Then 

(i) In the general settings TZs does not exceed Cph$~ d . 

(ii) In frames of proposition 4-11 ^3 does not exceed the right-hand expression of (4-78). 
(Hi) In frames of proposition 4-15(iii) H3 does not exceed the right-hand expression of 

(4.84). 
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(iv) In frames of proposition 4- 17(H) 7Z 3 does not exceed the right-hand expression of 
(4-88) as (4-75) is violated. 

(v) In frames of proposition 4-17(iii) 71% does not exceed the right-hand expression of 
(4-90) as (4-75) is violated. 

► This proposition together with corresponding results of section 4 implies all statement 
of theorems dealing with IZj rather than 1Z (theorems 0.9(i), 0.11(h)). It also implies 

Theorem 5.12. In frames of proposition 5.11 corresponding estimates hold also for Hi 
defined by (0.28), (5.61). 

It follows from proposition 5. 10(iii) that in frames of (0.33) with u(h) = h 2 we do 
not need any correction. On the other hand the following example shows the need in the 
correction term. 

Example 5.13. Consider the case q — 1, r = 2, 3 and fj = const and comeasurable and 
there are unremovable cubic terms in A so that the A\ = Ao + [i~ l M. with M. commuting 
with A (at each point x). Furthermore, let us assume that one of conditions is fulfilled: 

(5.62) 2 r = 2, f l= 2/ 2) 

(5.62) 3 r = 3, /i = / 2 + / 3 , / 2 // 3 . 

Then 

(5.63) 2 Ao = 2Z\Z X + Z* 2 Z 2 + D 2 5 , M = (ujZ\Z\) + (ujZ{Z\)\ 
(5.63) 3 Aq = hZ{Zx + f 2 Z*Z 2 + f 3 Z* 3 Z 3 + D 7 2 , M = {uZ;Z 2 Z 3 ) + (uZ\Z 2 Z z )* 

respectively with u = uj(x", x'", \i~ x h,D") 25 . 

We can already apply proposition 5.11, getting for cRi better estimate than estimate 
Cti2h l - d for n, but we can improve further the former one. Let us assume that \ou\ > e. 
Then (assuming large enough smoothness of symbols) one can apply a transformation by 
operator e th c with C = ^2-(3jZ*Zj with Hermitean (3j = f3j(x" \x"' ', [T l hD") and to 
transform operator (modulo terms containing factors fi~ l hDd or /i~ 2 ) to the same form but 
with Hermitean operator uj. So, M. = uAio with 



(5.64) Mi 



ZlZl + Z^* 2 r = 2 

Z{Z 2 Z ?J + ZiZ*Zl r = 3 



3 One should not forget to add Cji 2 Aq to avoid some non-semiboundedness related problems. 
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commuting with Aq. Then one can decompose with respect to common eigenf unctions of 
Ao and Aio and to break the systems into separate equations and one can easily prove that 

(5.65) III < C^~ d 

without any microhyperbolicity condition (which is much less than C^h 1 ~ d as ph <C 1). 
Let us select V = — Ylji^j + l)/?'/ 1 ^ an d H > h~z . 
Further, let us replace n(x, — z 2 , h) by 

(5.66) #{a : ^(2a 3 - + l)f 3 ph + V < 0} + m(x, h)6(p 2 - z 2 ) 

j 

where m(x, A, h) is an eigenvalue counting function for Ai restricted to the linear span 
of T Q such that J2j(^ a j + tyfj^h + V = 0. Then the error in will not exceed 

Cfiph 1 ~ d x (/ip 2 )" 2 . Taking p = pi we get this error not exceeding Cfx~^h~^~ d <C p^h l ~ d 
as long as ji hT*. 
But then (with 

(5.67) J £™™{x, 0)ip{x) dx = 

ffh r - d J j m(x, -e d , h)9(p 2 - £^(x) dxdi d + 0(Cfi-h-i- d ) 

which is exactly of magnitude pT* x p,h l ~ d = /i2/i 1_d and one cannot skip it without 
increasing the remainder estimate. 

5.9 Strong Magnetic Field 

Consider now case eh~ 1 \ \ogh\~ 1 < p < eh" 1 . The same arguments work again, with no 
exterior zone and necessity to join exterior and intermediate zones. Because of this one 
should deal directly with (5.5) rather than (5.24). 

Then (5.5) is preserved modulo 0(h s ) as T n is replaced by some lesser value which 
depending on assumptions is described in the previous subsections, then we apply method 
of successive approximations in the same way as above replacing U by U and ip by \psio 
and estimating an error in (5.5) cause by this; finally we replace T' n by T n again with 0(h s ) 
error and then T n by +oo. I leave details to the reader. 

Then one can easily prove 
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Proposition 5.14. Let eh' x \ log h\' x < p, < eh' 1 , (I, a) = (1,1) as q > 2, (I, a) = (1,2) 
as q = 1. T/ien one can rewrite (5.5) as J £ MW (x , 0)ip (x) dx with an error not exceeding 
Ch 1 " 11 asq>2, C^h 1 ' d asq = l 26 . 

► This statement combined with proposition 4.18 yield immediately theorems 0.13, 0.14. 

Proposition 5.15. Let eh" 1 ] \ogh\~ 1 < [i < eh" 1 , (I, a) >z (1,1) as q > 2, (I, a) h 

(1,2) as q = 1. T/ien under microhyperbolicity condition 0.1 one can rewrite (5.5) as 
J S MW (x, 0)i/j(x) dx with an error not exceeding Ch 1 ' d . 

► This statement combined with with proposition 4.19 yield immediately theorem 0.15. 

Proposition 5.16. Statements of propositions 5.8-5.11, theorem 5.12 and example 5.13 
remain true as eh~ 1 \ \ogh\~ 1 < p < eh' 1 . 

► This statement combined with propositions 4.20, 4.21 imply theorem 0.16. 
5.10 Superstrong Magnetic Field 

Consider now case eh' 1 < fi < C^h' 1 . The same arguments work again, with no interme- 
diate zone anymore. One can prove easily 

Proposition 5.17. Let eh' 1 <n< C^h' 1 , (/, a) = (1, 1) as q > 2, (I, a) = (1,2) as q = 1. 
Then one can rewrite (5.5) as J £ ww (x,0)ip(x) dx with an error not exceeding Ch 1 ' d as 
q > 2, Chh — d as q = 1. 

Proposition 5.18. Let eh' 1 < fi < C^h' 1 , (I, a) = (1, 1) as q > 2, (I, a) = (1,2) as q — 1. 
Then under microhyperbolicity condition 0.2 with I independent on r one can rewrite (5.5) 
as J £ MW (x, 0)ijj(x) dx with an error not exceeding Ch 1 ' d . 

Proposition 5.19. Statements of proposition 5.8 remain true as eh' 1 < \i < C^h' 1 . 

► These statements combined with proposition 4.22 imply theorem 0.16. 

26 One can improve it as there are no unremovable cubic terms but there is no point as the general 
remainder estimate 0{jih^~ d ) is larger 
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5.11 Ultrastrong Magnetic Field 

Consider now case ji > C^h' 1 . The same arguments work again, now applied to a single 
scalar equation: 

Proposition 5.20. Let \i > C^hr 1 , {I, a) = (1,1) as q > 2, (I, a) = (1,2) as q = 1. Then 
one can rewrite (5.5) as J £ MW (x,0)ip(x) dx with an error C jj r h r ~ d+1 as q > 2, Cjj r h r ~ d+ 2 
as q = 1. 

Proposition 5.21. Let \i > C$h~ l , (I, a) = (1,1) as q > 2, (I, a) = (1,2) as q — 1. T/ien 
under microhyperbolicity condition (0.56) one can rewrite (5.5) as J £ MW (x, O)ip(x) dx with 
an error not exceeding Cfi r h l+r ~ d . 

► These statements combined with proposition 4.23 imply theorem 0.17. 
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